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Trapped atomic ions have proven to be one of the most promising candidates for the re- 
alization of quantum computation due to their long trapping times, excellent coherence 
properties, and exquisite control of the internal atomic states. Integrating ions (quantum 
memory) with photons (distance link) offers a unique path to large-scale quantum com- 
putation and long-distance quantum communication. In this article, we present a detailed 
review of the experimental implementation of a heralded photon-mediated quantum gate 
between remote ions, and the employment of this gate to perform a teleportation pro- 
tocol between two ions separated by a distance of about one meter. 
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1. Introduction 

Quantum information research has the potential to drastically alter the fields of 
communication and computation. Efforts in quantum computation are driven by 
the prospect of using the features of quantum physics to tackle otherwise intractable 
computational problems. It was realized early on that controllable quantum systems 
may be used to simulate larger quantum systems far more efficiently than is possible 
using conventional computers P and that individual quantum systems might be used 
as quantum bits (qubits) for information processingP While both of these instances 
were crucial points, interest in quantum information increased dramatically in 1994 
when Peter Shor unveiled an algorithm that could be implemented on a quantum 
computer that would enable an exponential speed-up in the factorization of large 
numbersP Given that current encryption techniques, such as the RSA algorithm, 
rely on the relative inability of a conventional computer to factor large numbers, 
Shor's factorization algorithm linked quantum computation to the immediate and 
important issue of code-breaking. 

On the other hand, quantum information research has also yielded a new method 
of secure communication, where the security of the encrypted information is guar- 
anteed not by a lack of efficient mathematical or computational algorithms, but 
by the physical quantum properties of the information carriers. It has been shown 
that entangled pairs of qubits can be used to securely transfer information, where 
the presence of an eavesdropper is detected by the measurements made during the 
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protocoled The essential component for this protocol is the quantum no-cloning 
theorem, which proves quantum physics prohibits the cloning (copying) of an un- 
known quantum statePBy using a quantum communication protocol with pairs of 
correlated qubits, the back-action of measurement by an eavesdropper destroys the 
correlations, announcing their presence before any information is transmitted. 

Of course, the same features of quantum mechanics that are used to ensure the 
security of communication also prevent a simple "read-and-send" approach to trans- 
mitting quantum information. Any attempt to read or measure a single quantum 
superposition state results in a single measurement outcome, and thereby lacks the 
information needed to reconstruct the probabilities of the quantum superposition. 
While an estimation of an unknown quantum state could be obtained by simply 
measuring a large number of copies, the quantum no-cloning theorem forbids gener- 
ating identical copies of a single unknown quantum state. Nevertheless, a quantum 
state can still be transferred through the process of quantum teleportation.^ In 
the quantum teleportation protocol, a quantum state initially stored in system A 
can be recovered at system B without ever having traversed the space between the 
systems. The ability to teleport quantum information is an essential ingredient for 
long-distance quantum communication and may be a vital component to achieve the 
exponential processing speed-up promised by quantum computation. The essence of 
quantum teleportation lies in the non-local correlations, or entanglement, afforded 
by quantum physics. 

Quantum communication and quantum computation both use entanglement as 
an essential resource. Entanglement is the quantum correlations between systems 
that do not have well-defined individual properties. In addition to being useful for 
quantum communication and quantum computation, entanglement embodies the 
counter-intuitive depiction of nature predicted by quantum physics, and allows for 
explicit experimental tests of quantum theoryP Entanglement has been observed 
in a wide variety of systems, including: photons^ atomic ionsp^ superconduct- 
ing Josephson junctions)^ and neutral atoms in cavities,^ confined by an optical 
lattice,^ and in small ensembles.^ Hybrid systems composed of a single photonic 
qubit and a matter qubit (ionfSl ensemble) 15 l 16 l 17 IHl or atondES) have also shown 
entanglement. Here, we will focus on photons and atomic ions, both of which have 
already proven amenable to applications in quantum information. 

Atomic ions are one of the most promising systems for quantum information 
processing due to their long trapping times, excellent coherence properties, and 
exquisite control of the internal atomic states. To date, the largest entangled state 
of individually addressable qubits has been an 8-particle W-statepD realized in 
a system of trapped atomic ions that used the collective motion of the confined 
atoms to implement the entangling protocol. The experimental effort is now focused 
on scaling this system to larger numbers of qubits. One approach to scaling this 
systerrPHis to use microfabricated ion traps53Hl and advanced ion trap arrays H5EH 
Since these deterministic gate operations utilize the common modes of motion of 
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the trapped ions, an understanding of motional decoherence is also actively being 
pursuedHMMl 

Photons, on the other hand, are a natural choice of qubit for communication pur- 
poses, as they can quickly traverse the distance between locations with only small 
perturbations to the encoded quantum information. Already, a series of seminal ex- 
periments have used photons in quantum communication protocols over distances as 
large as 144 km .^Ultimately, though, the direct communication of quantum infor- 
mation over long-distances is impeded by the attenuation of light in air and optical 
fibers. For transmission through optical fibers, even at telecom wavelengths, which 
experience the least attenuation in fiber, the probability of transmitting a single 
photon over 1000 km is < 10~ 20 . Current fiber-optic information transfer mitigates 
the loss in signal amplitude by introducing repeaters along the transmission path 
to "boost" the signal along the way. Although application of the standard model 
of a repeater to quantum information is prohibited by the no-cloning theorem, an 
analogous "quantum repeater" has been suggested to enable the transfer of quan- 
tum information over arbitrary distancesl 31 l 32 lln this method, the distance between 
two points is broken up into a series of shorter segments, with a quantum memory 
at each connection point, or node. Entanglement can then be established between 
pairs of nodes, and subsequent segments connected via entanglement swapping, 
which is used to extend the entanglement over the entire length of the repeater. 
The final step is to use this long-distance entanglement as a resource to transfer 
quantum information over that distance by the process of quantum tclcportation. 

Integrating atomic ions (quantum memory) with photons (distance link) of- 
fers a unique path to large-scale quantum computation and long-distance quantum 
communication. Combining ions with photons enables long-distance quantum op- 
erations between stable quantum memories. These photon-mediated operations are 
mostly insensitive to the motional state of the ions, and thus can tolerate motional 
heating and do not require ground-state cooling. Moreover, the two-photon scheme 
presented here is not interferometrically sensitive to the optical path length dif- 
ference. Large-scale implementation may also be simplified by avoiding the need 
for complex trap arrays that allow for shuttling ions, or by enabling a higher-level 
architecture that links distant trap arrays. In addition, since these operations are 
mediated by photons, it may be possible to create hybrid matter quantum systems: 
interfacing atoms with solid-state qubits, such as quantum dots or NV centers. Fi- 
nally, while this photon-mediated operation is inherently probabilistic, it can still 
be efficiently scaled to enable generation of the large entangled states required for 
quantum information processing. 

In the following sections, we present a detailed review of the experimental im- 
plementation of a probabilistic quantum gate, and employment of this gate in a 
quantum teleportation protocol to transfer a qubit between two atomic ions sepa- 
rated by about one meter. We begin in Sec. [2] with the basics of trapping charged 
particles. Next we delve into the specifics of our atomic qubit, the ytterbium ion 
(Yb + ), in Seep] Section H presents the two-photon interference effect. In Sec. ^ 
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are the results of the experimental implementation of the heralded quantum gate. 
Finally, Sec. [6] illustrates the use of this gate in the long-distance teleportation pro- 
tocol. We conclude with an outlook of the potential practical application of these 
results to long-distance quantum communication and large-scale quantum compu- 
tation. 

2. Ion Trap 

The radiofrequency (rf) ion trap was invented by Wolfgang Paul, for which he 
shared the 1989 Nobel prizeP^ The first laser-cooling experiments with atoms were 
reported independently by Wineland et alJ^ using Mg + , and Neuhauser et alP^ 
using Ba + Q Since its introduction, the rf ion trap has been used for a wealth 
of applications, including atomic clocks, measurements of fundamental constants, 
mass and frequency spectroscopy) 36 ! 37 ! 38 ! 39 ! and quantum information science.^ 

According to Earnshaw's Theorem, "a charged particle cannot be held in a stable 
equilibrium by electrostatic forces alone. This result is also embodied in the 
freespace Maxwell equation V • E = 0, which in words states that electric field lines 
entering a region free of charges also need to exit the region; thus, any configuration 
of static fields will always be anti-trapping in some direction. Nevertheless, it is still 
possible to confine charged particles with electromagnetic fields by using either a 
combination of electric and magnetic fields (Penning trap) or dynamic electric fields 
(rf or Paul trap). In the following, we concentrate exclusively on the rf trap. 

We consider dynamic quadrupole potentials to trap charged particles. As an il- 
lustration in two dimensions, consider the hyperbolic electrode configuration shown 
in Fig. [l] which has the nearest electrode a distance R from the center of the struc- 
ture. Given that the left and right electrodes have potential Vbcos(f2:r*) applied, 
while the top and bottom electrodes are held at V, the potential between the 
electrodes is given by 

^ w = ^cos(Q r *)('l + * 2 ^y (1) 

which satisfies the boundary conditions, and Laplace's equation. The electric field 
produced by this potential is 

E(x,y,t) = ~V(j> h yp 

= -jg{xx - yy) cos(n T t) 

= -E (x,y)cos(tt T t), (2) 

where we have written the electric field as a product of its spatial variation Eq(x, y) 
and its time variation cos (f2yt). The force in the ^-direction on a particle of mass 

a The experiment by Wineland et al. used a Penning trap to confine the ions. Subsequent experi- 
ments by Wineland et al. have used the type of rf trap presented here. 
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Fig. 1. Hyperbolic electrodes. This two-dimensional view shows hyperbolic electrodes where the 
green electrodes are held at ground, and the red electrodes are at Vb cos(fi;ri). The resulting 
potential is shown as a contour plot between the electrodes (shades of blue) for t = 0. 



m and charge e is then 



F x = mx = — ^-^ cos(flTt)x (3) 
R 2 



yielding the equation of motion 

A + As? cos(fl T t)x = 0. (4) 

This is actually just a simplified version of a Mathieu equation. The general form 
of the Mathieu equation is: 

d 2 u 

— + (a u + 2q u cos(2r)) u = 0. (5) 
dr 2 

Equation [4] can be put in this form by making the substitution 2r — ilxt. By the 
chain rule, we then have = £p ^ = ^ £p , so that Eq. [4] can be written as 

d 2 x 

j^+2q x cos(2T)x = 0, (6) 

where in the last line we have defined q x = (2eVo)/(mR 2 ^l'^,). Note that compared 
to the general form of Eq. [5| here we have a x — 00 

The Floquet Theorem suggests a general solution of the forrrP^l 
00 00 
x(t) = A J2 C 2n cos((2n + f3)T) + iB C 2n sin((2n + (3)t) 7 (7) 



n— — 00 n— — 00 



b In the present case, we have neglected the third dimension of the problem. The ion can be 
confined in the third dimension using static potentials (see below), in which case a x ^ 0; a good 
review is Ref. 41 
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where A, B are determined by the initial conditions. We can solve for (3 and the 
coefficients C 2n by plugging Eq. [7] back into Eq. [6] We then find 

oo 

A Y C 2 „(2n + /3) 2 cos((2n + /?)r) 

71 — — OO 

OO 

+iB Y C 2n (2n + /3) 2 sin((2n + /3)r) 

n— — oc 

OO 

= 2q x A Y C 2ri cos(2r) cos((2n + (3)t) 

71 — — - OO 
OO 

+i2q x B Y C , 2 „cos(2r)sin((2n + /3)T). (8) 

71 — — OO 

Using the trigonometric product-to-sum relation for cosing] we get 

OO 

A Y C 2n (2n + /3) 2 cos((2n + (3)t) 

71= — OO 
OO 

+iB Y C 2n (2n + P) 2 sm((2n + f3)r) 

71 — — OO 

OO 

= ^ C 2 „ (cos((2n + /3)r + 2r) + cos((2n + /3)r - 2t)) 

n— — oo 

OO 

+i2q x B Y C 2 „cos(2T)sin((2ri + /3)r) 

71 — — OO 

OO 

= q x A Y C ^n-2 cos((2n + (3)t) 

71— — OO 
OO 

^ C 2n+2 cos((2n + /3)r) 

71 — — OO 

OO 

+i2q x B Y C , 2 „cos(2r)sin((2n + /3)r), (9) 



71 — — OO 



where in the last step we just altered the indexing (possible since the sum goes 
between ±oo). This allows us to easily match the cosine terms, yielding a recursion 
relation: 

(2n + 8) 2 

- K 2n C 2n + C 2n _ 2 + C 2n+2 = 0, with K 2n = 1 ' ' ' . (10) 

q x 

This recursion relation allows us to calculate some useful relations. By setting n = 

c The product-to-sum relation is cos(cr) cos(/3) = \ (cos(a + 0) + cos(q — /3)). 
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in Eq. 10 we obtain an expression for j3: 



Qx L^O 



Rearranging Eq. [TUJ we also find 



C 2 „+2 K 2n - % 



(12) 



c 2 „ 

and we can then plug Eq. [12] back into itself recursively, and obtain the continued 
fraction expression 

C2n 1 



C211+2 K 2n 



(13) 



K 2n - 2 - J - 



Similarly, we also get 



f ' 2 " ' , • (14) 



C2n-2 Kl n — 



L 2n — ~T? r 

"271 + 2 — — 

Plugging Eqs. [14| and |13"| into Eq. [TT] allow the calculation of (5 to any order in q x . 
Explicitly, we have 

P 2 = q x 



C_2 C*2 

, Co C 
1 



K« - — J-t- X 2 - 



(15) 



if_4-— 2 K4-- 



The solution to the Mathieu equation has been shown to be stable for q x less 
than about 0.9P^To lowest order in q x , we find 

' 1 1 



and thus 



if-2 K 2 

Mt + tJ 

f, (16) 

Given the above expressions for /3 and the coefficients C2n, we can solve for the 
trajectory of the ion, x(r). Assuming ^ < 1, we will take C±4 ~ 0. By assuming 
the initial condition B — 0, from Eq. [7] we then obtain 

x[f) w AC cos(/3r) + AC_ 2 cos((/3 - 2)r) + ^4C 2 cos (03 + 2)r) . (18) 

Since g x -C 1, Eqs. |13| and |14| can be approximated as 

^ C+2 C+2 C+2 4C+2 
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where in 
09 T 2) 2 : 



x(t) 



the second-to-last step, we used the approximation of Eq. 17 and 
s 4. Plugging this value for C±2 into Eq. 18 yields 

Qx 
4 

2 

2 



AC cos(/3r) 
AC cos(/3r) 
AC cos (/3t) 



(cos((/3 - 2)r) + cos((/3 + 2)r)) 
cos(/3t) cos(2r 
■ cos(2r) 

We can finally put this back in terms of t by recalling our earlier substitutions: 

Qx 



1 



(20) 



x(t) « ACo cos 



AC cos 



= ACq cos(uj x t) 



q x 

2 1 / 2 mn T R 2 

q x 



1 



cos(riri) 



1 



<7.r 



cos(f2j^) 



1 



COs(fiTi) 



In the last step we have defined the "secular frequency" of the ion as: 

eV 



(21) 



(22) 



x 2V2 m o Tj R2 • 

The secular frequency characterizes the relatively slow harmonic motion of a 
charged particle confined by an oscillating quadrupole field. In analogy, the har- 
monic motion executed by the charged particle can be derived from an effective 
harmonic pseudopotential produced by the quadrupole fieldP^ 



(23) 



The secular motion of a collection of ions has been utilized to perform deterministic 
quantum information processing tasks in several experiments.^ 

Superimposed on the secular motion in Eq. [21] is a modulation at the driving 
frequency of the trap, f^T- The oscillation at fi-r is known as the "micromotion," and 
has an amplitude proportional to the distance from the center of the quadrupole 
field. Given that q x <C 1, the amplitude of the micromotion is typically much smaller 
than the secular motion. However, if a static offset field causes the secular motion 
to be centered a considerable distance from the center of quadrupole field, then it 
is possible for the amplitude of the micromotion to be much larger than the secular 
amplitude. In practice, static offset fields are carefully compensated to ensure the 
particle executes its secular motion about the center of the quadrupole fieldP^ 

A representative picture of the motion of a particle confined by the two- 
dimensional quadrupole field of Eq. [2] is shown in Figure [2] In this case, we have 
taken R = 0.46 mm, n T /{2n) = 38 MHz, m = 171 amu, and V = 1 kV. By the 
equations above, we see this results in u} x /(2tt) — 1.3 MHz and q x = 0.1. 

The derivation above can be extended to confinement of the charged particle 
in three dimensions by either the addition of another rf electrode along the z-axis 
or by capping the third dimension with static potentials. In the latter option, the 
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Fig. 2. The motion of a trapped charged particle, as given by Eq. |21| with R = 0.46 mm, 
Oy/(27r) = 38 MHz, m = 171 amu, and Vo = 1 kV. The slower, larger amplitude oscillation at 
uj x /(2tt) = 1.3 MHz is the secular motion of the particle, resulting from the time-averaged force 
of the inhomogenous electric field (red dashed line). Superimposed on the secular motion is the 
micromotion that occurs at the drive frequency of trap (Qt I (2tt) = 38 MHz) with amplitude 
proportional to the excursion of the particle from the center of the quadrupole field (blue solid 
line). 



hyperbolic electrodes can be extended along the z-axis to produce an rf nodal line 
in this dimension. By making the static confinement of the particles along z weaker 
than the rf confinement in the x, y directions, the charged particles can be arranged 
in a linear crystal. In such a "linear trap" all of the particles execute their secular 
motion about the center of the quadrupole field. The addition of a static potential 
along z produces a potential near the center of the trap of the form)^ 

Uo ( 2 x 2 + y 2 



i\ z 2 



2e 



* 2 -^!]. (24) 



Here we have defined the distance to the static electrode as zq, and the oscillation 
frequency along z as lu z — ^/ 2eXJ q j (vnzfy . In addition to confining the particle 
along z, this static potential also alters the potential in x and y, so that the secular 
frequency of the ion becomes 



(25) 

Thus, the static potential along z weakens the tranverse (x,y) confinement of a 
charge particle. However, the rf potential is typically much larger than the applied 
static potential, so the trap remains stable. 

The two-dimensional hyperbolic electrode structure presented in the previous 
section produces the ideal quadrupole potential for trapping charged particles. Of 
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course, in practice we are often required to alter the structure of the trap electrodes 
from the ideal hyperbolic structure to conform to other design parameters, such as 
optical access, multiple trapping zones, or simplicity of construction. Nevertheless, 
as long as we retain the overall symmetry presented above, then near the trap center 
the potential can be approximated by the hyperbolic potential given in Eq. |lp^ 
The effect of non-hyperbolic electrodes can be characterized by the addition of a 
geometric scaling factor rj sc (generally of order unity) , so that the modified potential 
near the center of the trap is approximately represented by!^ 



Since rj sc is just a constant factor, it carries through the rest of the equations derived 
in the previous section. 

Numerical simulations of the candidate ion trap can be used to ensure the ex- 
pected properties are consistent with the requirements of the experiment and the 
capabilities of the available equipment. Given that the wavelength of oscillation of 
the potential applied to the electrodes (order of meters) is typically much larger 
than the size of the trap (< 1 mm), electrostatic simulations are sufficient. As an 
example, consider the four-rod trap illustrated in Fig. [3] this trap design is used in 
all of the experiments presented in the following sections. The simulated quadrupole 
potential for the four-rod trap is shown in Fig. |4ja) , with the resulting pseudopoten- 
tial illustrated in Fig.[4jb). Comparing Fig.[3Ja) with the ideal quadrupole potential 
of Fig. [I] we see that near the center of the trap the four rods produce a potential 
nearly identical to that generated by hyperbolic electrodes. 



Fig. 3. Four rod trap. The trap design shown here is the type used for all experiments presented in 
this work. Two of the four rods have radiofrequency (rf) potentials applied to create an oscillating 
quadrupole field for confinement in the x,j/-plane. The other two rods are held at ground. Static 
potentials are applied to the needle electrodes to provide confinement along the 2-axis of the trap. 

The motion of a charged particle in a trap can be viewed as a three-dimensional, 
uncoupled harmonic oscillator along the principal axes of the trapping potential. 
In other words, motion of the particle along one of the principal axes of a trap is 
independent of the other two principal axes. Therefore, the principal axes of the 
trap present a natural coordinate system for the electrode structure. In addition, 
knowledge of the orientation of the principal axes is vital for efficient laser cooling; 
if the incident light is perpendicular to one of the principal axes, the particle will 




(26) 
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Fig. 4. Numerical simulation of four-rod trap potentials, (a) The quadrupole potential produced 
by application of voltage to the two rf rods, while the other electrodes are held at ground (shown 
for the point of maximum amplitude in the oscillation), (b) The pseudopotential derived from 
the quadrupole potential of (a) , using Eq. |23| The numerical simulation was performed using the 
boundary element method (BEM) electrostatic modeling software CPO-3D from Charged Particle 
Optics, Ltd. 



not be cooled along that direction. While in simple trap structures such as our 
four-rod trap, the orientation of the principal axes is clear from the symmetry of 
the electrodes, in general we determine the principal axes of the trap by use of the 
Hessian matrix of the potential, which in two-dimensions isP^ 



(27) 



Here <fi S im is the simulated trap potential (pseudopotential plus static), and xg and 
yo are the coordinates where the Hessian matrix is evaluated, which in our case 
would be the center of the trapj^] The eigenvalues of this matrix are related to the 
angle by which the principal axes are rotated with respect to the coordinate axes 
(used in the partial derivatives). The Hessian matrix seeks out the directions of 
greatest and least curvature, which are precisely the principal axesQ 

The four-rod trap can now be evaluated for specific parameters, and compared 
to the ideal hyperbolic case to determine the geometric scale factor, r/ sc . Taking 
V rf = 1 kV, V dc = 80 V, m = 171 amu, R = 0.5^2-0.25 = 0.46 mm, and fl T = 38 
MHz, a numcrcial simulation determines a secular frequency uj XyS i m /{2'K) = 1.17 
MHz for the four-rod trap. Plugging these same parameters into Eq. [22] for the 
ideal hyperbolic trap, we find w x ,hyp/{^) — 1-26 MHz. Thus, the geometric scale 

d In a linear trap, the third principal axis is always clear from structure of the trap. In Fig. [3] it 
is parallel to the rods of the trap (along the line defined by the needle electrodes). 
c Degeneracy of the principle axes can be broken by the addition of a static potential along either 
x or y, as noted in Ref. 1411 
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factor is r) sc = oJ x ,sirn/^x.hyp = 0.93. The total depth of the trap can also be 
determined by the simulation from the depth of the pseudopotential-well, and for 
these parameters we find the trap depth to be approximately 10 eV. This trap 
depth corresponds to a temperture of more than 10 5 K, reiterating the fact that 
the rf trap provides excellent confinement of charged particles. 

3. Ytterbium Ions 

Trapped atomic ions have long been recognized as a promising implementation 
of quantum bits (qubits) for quantum information processing P3Hdue in part to 
long trapping lifetimes, long coherence times of particular internal electronic states, 
and the exquisite control attained over both the internal and external degress of 
freedom. The hydrogen-like ions that have been directly cooled and manipulated for 
applications in qu antum informa tion i nclude Ba+) 48 l 49 l 50 l Be+J^ 1 Ca+ j52|53|54|55|56j 
Cd+J^Mg+JHHSr+J^and Yb+) 60 l 61 lln Table [l] various properties of some atomic 
ions are compared. 



Table 1. Ion qubit comparison. The NIST database^] was used to determine many of the numbers quoted here. 





Be+ 


Mg+ 


Ca+ 


Zn+ 


Sr+ 


Cd+ 


Ba+ 


Yb+ 


Hg+ 


isotope 


amu) 


9 


25 


40, 43 


67 


87, 88 


111, 113 


135, 137, 138 


171, 173 


199, 201 


nuclear 


spin 


3/2 


5/2 


7/2 


5/2 


9/2, - 


1/2, 1/2 


3/2, 3/2, - 


1/2, 5/2 


1/2, 3/2 


2 S 1/2 hfs (GHz) 


1.25 


1.8 


3.2 


7.2 


5, - 


14.5, 15.3 


7.2, 8, - 


12.6, 10.5 


40.5, 30 


P is (THz) 


0.2 


2.75 


6.7 


26.2 


24 


75 


50.7 


100 


274 


2 S 1/2 <-> 


2 P , 


313.2 


280.4 


397 


206.3 


421.7 


226.5 


493.5 


369.5 


194.2 


2 s 1/2 ~ 


2 Pz/2 


313.1 


279.6 


393.5 


202.6 


407.9 


214.5 


455.5 


329 


165 


2 £>3/2 <-> 


2 Pl/2 






866.5 




1091.8 




649.9 


2438 


10747 


2 £>3/2 *-* 


2 P 3 /2 






850 




1003.9 




585.5 


1350 


994.7 


2 D s/2 ~ 3 


[3/2] 1/2 
















935.2 




2 D 5/2 <-> 


2 P , 
P3/2 






854.4 




1033 




614.3 


1650 


398.5 



The ytterbium ion (Yb + ) has several advantages. The strong 2 Si/ 2 <-> 2 P\/2 
electronic transition near 369.5 nm is suitable for use with optical fibers, making 
schemes that require the coupling of atomic (hypcrfine) qubits to photonic (opti- 
cal) qubits feasible.^ Moreover, the large fine structure splitting of Yb + makes it 
amenable to fast manipulation with broadband laser pulses EMMEZI Finally, the 
spin-1/2 nucleus of 171 Yb + allows for simple, fast, and efficient preparation and 
detection of the ground state hyperfine levels.^ 

In the experiments presented in the following sections, an ytterbium ion (Yb + ) 
is produced by photoionization of neutral ytterbium using a two-photon, dichroic, 
resonantly-assisted process. The resulting Yb + atom is confined in the four-rod 
rf trap described in the previous section. The ion is Doppler-cooled by light at 
369.5 nm, which is slightly red-detuned of the 2 S*i/2 <-> 2 -Pi/2 transition depicted 
in Fig. [5] Efficient cooling of the atom requires addtional light, because the 2 Pi/2 
state also decays to the metastable 2 £*3/2 level with a measured probability of about 
0.005P Light at 935.2 nm is used to drive the atom from the 2 D 3/2 to the 3 [3/2] 1/2 



Quantum Logic Between Distant Trapped Ions 13 



level, from which it quickly returns to the 2 6*1/2 ground stated An additional 
complication arises from the presence of the low-lying 2 F 7 / 2 state. Despite the fact 
that there are no allowed decays from the four levels used in cooling to 2 F 7 / 2 , 
the ion falls into this state a few times per hour, probably due to collisions with 
residual background gasl 69 l 70 HU Laser light near 638.6 nm depopulates the 2 F 7 / 2 
level, returning the ion to the four- level cooling scheme. Trapping times of several 
months have been observed. 



^ 1 [5/2] 5 , 



N 

I" 



TT 3 [3/2] 1/2 



■'312 
T - 52.7 ms 



CD 
00 
CD 



T= 7.2 ms v??, 



' 7/2 
T- 5.4 years 



Fig. 5. Relevant levels of the Yb + atom (to scale). Transitions shown with solid lines are driven by 
laser sources in the experiment. Numbers given in parentheses are branching ratios. The lifetimes 
of some of the excited states are also given. Measurement of the lifetime and branching ratio for 
the 2 Pi/2 level is given in Ref. 1 72 1 and Ref. 1611 respectively. Wavelengths for decays shown in gray 
are from Ref.|62] Lifetimes of the 2 D 3/2 , 2 D 5/2 , 3 [3/2] 1/2 , and 2 F 7/2 levels are from Ref s. 17311741 
1751 and 1761 respectively. Branching ratios out of 3 [3/2]jy 2 and 2 ^3/2 are from Ref. 1771 while the 
2 D 5 / 2 branching ratio is from Ref. 1741 



An external magnetic field of about 5 gauss is applied to provide a quantization 
axis for definition of the polarization of the photons emitted by the atom, break 
the degeneracy of the atomic states, and suppress coherent dark state trapping.^ 
The qubit states are chosen to be the first-order magnetic field-insensitive hypcrfinc 
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"clock" states of the 2 S l 1 / 2 level, \F = 0, rap = 0) and \F — l,mp — 0), which are 
separated by 12.6 GHz and defined to be |0) and |1), respectively; here, F is the total 
angular momentum of the atom, and wip is its projection along the quantization 
axis. Applying light at 369.5 nm resonant with the 2 S±/2\F = 1) «-> 2 P\/2\F — 1) 
transition initializes the qubit to the state |0) by optical pumping. Subsequently, any 
combination of |0) and |1) can be prepared by application of microwave radiation 
at 12.6 GHz with controlled duration and phase. The state of the atomic qubit 
is determined by using standard fluorescence techniques. By illuminating the ion 
with light at 369.5 nm resonant with the 2 Si/ 2 \F = 1) «-> 2 P\/2\F — 0) transition, 
an ion measured to be in the state |1) scatters many photons, while for an ion in 
the state |0) the light is off-resonance and almost no photons are scattered. We can 
then discriminate between the two qubit states with better than 98% fidelity.^ 



4. Two Photon Interference 

The interference of identical single photons was first observed by Hong, Ou, and 
Mandelf^ and Shih and AlleyP^ In these experiments, pairs of photons created by 
parametric down-conversion were directed to interfere at a beamsplitter. It was ob- 
served that the interference of these identical photons at the beamsplitter resulted 
in two simultaneously impinging photons always exiting the beamsplitter by the 
same port. Thus, coincident detections behind the beamsplitter were highly sup- 
pressed. This quantum interference effect is an essential component of the remote 
atom entanglement protocols reviewed in the subsequent section. 



4.1. Theoretical interference of photons at a beamsplitter 

We consider photons impinging on a 50:50, nonpolarizing beamsplitter, such as the 
one shown schematically in Fig. [6] We can describe photons and the effect of the 
beamsplitter using the usual photon creation operator a' and the vacuum state |0). 
A single photon in mode n (as in Fig. [6]) can then be written as |1„) = aJjO). The 
operation of the beamsplitter is then encompassed by^D 

4 = ^|(4 ~4) 

4 = ^(4 + 4)- (28) 

Of note is the sign difference between the decomposition of modes 3 and 4 in terms 
of modes 1 and 2. The minus sign is the result of the 7r-phase shift experienced 
by light reflected at an interface where the index of refraction changes from low to 



high 82 ' and ensures energy conservation. Equation 28 may instead be expressed as 

4 = -^(4 + 4) 

4 = -J= (-4 + 4), (29) 
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beamsplitter (BS) 




Fig. 6. Schematic of a 50:50, nonpolarizing beamsplitter (BS). Modes 1 and 2 are the two input 
paths to the beamsplitter, while modes 3 and 4 are the exit paths from the beamsplitter. 



which is effectively a decomposition of incoming photons in terms of the exiting 
photons. 

The above formalism makes it easy to see the effect of the beamsplitter on 
impinging photons. We first consider the trivial case of a single photon in mode 1. 



lit) = 410} 
i 

= -L(|i 3 o 4 ) + |o 3 i4». 



(4 



4)|o) 



(30) 



As expected, there is an equal (50%) probability of detecting the photon in mode 
3 as in mode 40 With only a single photon impinging on the beamsplitter, there is 
zero probability of getting a detection in both mode 3 and mode 4 simultaneously. 

We next consider two photons impinging on the beamsplitter: one in mode 1 
and one in mode 2. Following the same procedure as above, we then find 



lila) = 



it At I 



(4 + a*) (-4 + 4)|o) 



= *((4) 2 -(4) 2 )|o) 



'On a sidenote, notice that Eq. |30| appears to indicate that the effect of the beamsplitter is to 
produce an entangled state, yet only one photon is involved. The remedy to this apparent paradox 
is to realize that the entanglement here is betwee n two modes (mode 3 and mode 4) p3 although 
there is some controversy over this interpretationEHU] 
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= -^(|0 3 2 4 )-|2 3 04» (31) 

and again there are no simultaneous detection of photons in modes 3 and 4. Instead, 
both photons either go into mode 3, or both photons go into mode 4, with equal 
probability. This simple analysis qualitatively explains the two-photon interference 
effect observed by Hong, Ou, and MandelJ^ and Shih and AlleyP^ 

Now we generalize this formalism in a way that enables us to describe the ideal 
expected signal of the experiment. We closely follow the treatment presented in 
Ref. IHH but expand it to explicitly derive the signal expected from our experimen- 
tal implementation. As described in the following subsection, the experiment to 
demonstrate this two-photon interference effect will consist of repetitive fast ex- 
citation of two trapped atomic ions, and detection of the spontaneously emitted 
photons behind a beamsplitter. If we assume that the repetitions of the experiment 
are well-separated (as, indeed, is the design), then we can write the electric field 
operator of the mode j for the rfi 1 repetition as 

Ef )+ (t) = ^(r^e-i^^eit-nt^af, (32) 

where t p is the time between repetitions, r is the natural lifetime of the excited 
state, and is the Heaviside step function. The exponential decay factor accounts 
for the probability to detect a photon a given time interval following excitation. 
The 1/-\A" factor ensures that the total probability of detecting a photon is unitless 
(integral over time of Eq. [35] below). The factor A is an amplitude that depends on 
the photon generation and collection efficiency. The function £ (r, t) describes the 
spatial mode of the photons, and can also be used to account for temporal offsets. 
For the remainder of this section, we ignore both A and £(r,i)[s]The full electric 
field operator for the mode j is then simply 

^(t) = X>j n)+ W- (33) 

n 

The electric field operator describes the probability of detecting a photon at any 
given time. 

The A^+l photons produced in the mode j by N+l repetitions of the experiment 
can be represented by the wavefunction 

(n a j n)t ) ( 34 ) 

Vn=0 / 

where we have assumed that repetitions are sufficiently separated in time so that 
overlap of photons produced by different repetitions of the experiment is negligi- 
ble; this allows us to write photons produced by different repetitions as separate 
"repetition modes" denoted by n with the same spatial mode j. 



g We will find these factors to be useful later, when we calculate the error in the fidelity of the 
quantum gate and teleportation protocol due to mismatch of the photon spatial modes at the 
beamsplitter. 
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If a train of photons is produced in mode j, then the probability of detecting a 
photon in mode j at time t is simply given by the first-order correlation function 

^MErm+m^ = <0| (f[ af) (j2 j^e-^- kt ^&(t-kt p )afA 

\n=0 / \fc=0 ^ T ) 
x (E ^e'^-^y^{t-mt p )a\ m) \ (f[&?A |0) 

\m=0 ^ T ) \n=0 / 

1 N 

= - e-5(*-fe* P )/^ e -5(*-™* P )/^e(t- kt p )e(t-mt p ) 

k,m— 

K<oifn^)as fc,t ^ m) fnaS B)t ) i°) 

\n=0 / \n=0 / 

1 N 

= - V e-3(*- fc *^/ r e^2(*-« l *p)/-e(i-A:t p )e(t-mt p )4r„ 

T ^ — ' 

k,m— 

1 w 

= - e-^ kt "^ T e(t - kt p ), (35) 



fc=0 



where, as mentioned above, we have ignored both A and £(r, t) from Eq. 32 Figure|7] 
illustrates this result for the case N = 10. As shown in the figure, the probability of 
detecting a photon from the wfi 1 repetition of the experiment decays exponentially 
with decay constant given by the natural lifetime of the excited state (t). The 
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Fig. 7. Theoretical single photon first-order correlation function, given by Eq. |35| This shows the 
relative probability of detecting a photon in mode j at time t following a series of 11 excitations 
(N = 10). The charateristic exponential decay after each repetition/excitation in the experiment 
is determined by the natural lifetime of the excited state of the atom producing the spontaneously 
emitted photons. The repetition rate shown is 20 times the natural lifetime (r). 
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probability for a joint detection in both mode 3 and mode 4, given a train of photons 
being produced in mode 1, is given by the second-order correlation function 

Pj{t,t d ) = 1 {^ N \E^{t)E^{t + t d )E+{t + t d )E+{t)\xl; N ) l 

= (0| ( ft a["A (f2^e-^-^e(t-kt p )a^A 

\m=0 J \k=0 V / 

x (j2 -^e-^ t+t ->- nt *y T G(t + t d - nt p )atA 

\n=0 ^ T ) 

x (E ^^^-«/ T 8(t + 1„ - r^a^j 

x IE -^e-i(*-^0(t - *tp)4« ) ( ft 4 m)t ) |0>, (36) 

\fc=0 ^ T ) \m=0 / 

where ^ is the time delay between the two detections. The creation and annhilation 
operators for modes 3 and 4 can be decomposed in terms of modes 1 and 2, using 
Eq. [28] Since in the above equation there are no photons in mode 2, it is clear that 
any term with an annihilation operator of mode 2 will immediately vanish. This 
allows us to simpifly Eq. [36] to obtain 

p j(t,td) = -L<0| (n a^A (j2e-^-^e(t-kt p )a[ k A 
\m=0 / \k=0 I 

x (j2 e-*< t+t "-' rt *>/ T e(f + < d - ni p )a< n)t J 

\n=0 / 

x ( ^ e-iCt+td-ntpJ/rg^ + <rf _ ntp ) a W j 

\n=0 / 

x ^e-i^-^eCt-fctpJSi* 5 J ( f[ 4 ro)t ] |0) 

Vfc=0 / \m=0 / 

= ^ E E e-^-^^^e-^+^-^^^G^ - Wj,)0(t + t d - Trip) 

k=0 n=Q 



n»(m) A »(fc)t-.(n)t~(n)-.(/s) | TT »(m)t \ i nX 
m=0 / \m=0 / 



iV iV 



^ ^ ^ e -(t-fctp)/r e -(t+t d -r ltp )/r e ^ _ ItfJQQ. + frf _ ^ 



4r 2 

fc=0 n=0 



x (1 - 4„) , (37) 

where the last step is simply a consequence of having an annihilation operator for 
both k and n, but only one creation operator for m; thus, if k = n the term goes to 
zero. Of course, we are interested in the total probability of a joint detection as a 
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function of the time delay td between the two photon detections. In order to obtain 
this function, we integrate Eq. 37 over t: 

/oo 
Pjoint(t>td) dt 
-oo 

N N 

/oo 
e -(t-kt p )/T e -(t+t d -nt p )/T Q ( t _ ktp }Q( t + U _ ntp j dt 
-oo 
1 N N 

= 5 \EE(i-4„) e - (td - t " (fc+ " ))/T 

k=0 n=0 



x e~ 2t/T 0(t-kt p )0(t + t d -nt p ) dt 



OO 

N N 



1 -2T,J2( 1 ~ s ^ e ~ {td " tp(k+n))/7 



4 T 2 

fc=0 n=0 



oo 

x / e~ 2t/r e(£ + t d - nt p ) dt 



N N 



4r 2 

e - 2u/T Q(u) du. (38) 

td+t p (k— n) 

In the last line we've set u = t + td — nt p . The integration is completed by looking 
at two cases: 



P.jtot{td < -t P {k - «)) - A E E (! - <W e(*«« +t '( fc - n »/' 



x / e- 2u / T Q(u) du 

N N 







^^(l-4n)e^ +f '« 



fc=0 n=0 

X 

AT N 



AT AT 



fc=0 n=0 



and 



N N 



fe=0 n=0 
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I e - 2 "/ T 6(u) du 

i+t p (k~n) 
N N 



f 

1, . 

^Y,Y,( 1 - s ^ eitd+tp{k ~ n 

i 

-EE( 1 -^) e " (td+tp(fe " n))/r 



k=0 n=0 

\-L e -*»/r 

2 \ td +t p {k~n) 
N N 



k=0 n=0 

By inspection, Eg. [39] and [40] can be written succinctly as 

N N 



(40) 



(41) 



fc=0 n=0 



This result is plotted in Fig. [8] As expected from Eq. 30 the probability of simul 
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Fig. 8. Theoretical single photon second-order correlation function, given by Eq. |41| This shows 
the relative probability of detecting a photon in both mode 3 and mode 4 as a function of the 
time delay (t^) between the two photons for N = 4 (red, dashed) and N = 20 (blue, solid). The 
repetition rate shown is 20 times the natural lifetime of the atomic excited state (r). The decrease 
in the peak amplitude of the probability of joint detection for increasing delay time results from 
having a finite train of impinging photons. The lack of a peak at delay time = is characteristic 
of a single-photon source. 



tancously detecting a photon in both mode 3 and mode 4 is zero, indicated by the 
lack of a peak at delay time td — 0. Experimentally, the number of detection events 
at td = relative to the height of the adjacent peaks can be used to evaluate the 
fidelity of single-photon generation ) 86 ! 87 ! 

We now look at the signal we would expect from photons impinging on both 
input ports of the beamsplitter. Equation |34| is used to express the photonic wave- 
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function for a train of pairs of photons in modes 1 and 2 as 

\*N)i*=(f[ti n A(fi£ m A io) 



/ \m— 

/AT \ 

4 (n)t-(n)t 



= ^n«r«rjio>, (42) 

where the second line follows from the fact that [a^^a^" 1 ^] = V m,n. Here we 
assume that the "repetition modes" are common to the two pulse trains; in other 
words, we assume to have pairs of identical photons impinging on the beamsplitter]^] 
The second-order correlation function for pairs of photons can then be calculated 



in a way similar to Eq. 36 



P2 P h,j(t,t d ) = i,2^ N \E^{t)E^{t + t d )E+{t + t d )E+{t)\^ N ) ia 

= (o| ( ft ai m) 4 m) ) (j2^e-^-^e(t-kt p )a^A 

\m=0 / \fc=0 ^ T ) 

\n=o v r / 

x ( £ -^e-^W-^-ntp)/^^ + f d - nt p )4" } ) 
\n=0 v T / 

\k=0 v T / 

xfn^ m)t 4 m)t )|0). (43) 

\m=0 / 

Of course, the time dependence of Eq. [43] appears identical to Eq. [37] Thus, the 
integration over time will be identical to the single-photon case, and so for the 
present calculation we only need to focus on the effect of the additional creation 
and annihilation operators. Writing the operators for modes 3 and 4 in terms of 
modes 1 and 2 yields 

/ X (mWmA ~(k)U(n)U(n)~(k) ( A A (m)t^(ro)t | | n \ 



iW(Ml) oc (0| J] a i 4 a^a^a^a^ JJ a[ 

\m=0 / \m=0 

1 /n i ( TT ~(m)-(m)\ /»(fc)t ~(fc)t\ /-Wt , 



\m=0 



h Slight temporal mismatch between photons in modes 1 and 2 could be accounted for in the 
definition of the electric field operator, Eq. |32| Since at the present we are interested in deriving 
the ideal signal that would be measured by this kind of photon interference experiment, we assume 
the ideal case of photons from the same repetition having perfect temporal overlap. 
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\m=0 / 
1 /nl I TT -(m)-(m) 1 /-(fc)t-(n)t-(n)-(, 

-(0| 1J_ a\ >a y 2 \a\> a\ ' a\ ' a\ 



(k) ~(k)U(n)U(n)~{k) 

— a a) a, a) 



\m—Q 



. -(fc)t~(n)t~(n)~(fc) . »(fc)t-.(n)t-(n)-(fc) » (fc)t - (n)t - («) - O) 

+a 2 ' a-[ a]^ a 2 + a\ a 2 y a 2 'a^ — a\ a 2 a 2 



+a^a^a^aP) ]J a^a™ |0> 



Vm-0 



= (1 - 4n) - fan + 1 + 1 - fan + (1 - <5fcn) 

= 4(l-5 fc n), (44) 

where the sum over A: and n is implicit. The final solution, after integrating over t 



(with integrals identical to Eq. 38 1, is 



N N 

P2 P h,Jtot(t d ) ^EEt 1 - e-l^ +t ^ fe -")l/ T , (45) 



2r ■ 

fe=0 n=0 



which is just Eq. 41 multiplied by a factor of 4. This indicates that there is zero 
probability of obtaining simultaneous dcctections in both modes 3 and 4 if the 
incoming photons in modes 1 and 2 arc identical, as was determined qualitatively 



in Eq. 31 The additional factor of 4 results because in this case one photon in mode 
3 and one photon in mode 4 separated by non-zero delay time td can be the result 
of: 2 sequential photons in mode 1; 2 sequential photons in mode 2; 1 photon in 
mode 1 followed by 1 photon in mode 2; 1 photon in mode 2 followed by 1 photon 
in mode 1. 

Note that the fidelity of a two-photon interference experiment can be evaluated 
by taking the ratio of the number of detection events at t$ = to 1/2 the number 
detected at a multiple of the experiment repetition rate td = nt p . The factor of 
1/2 is derived from the second-order correlation function for an incoming train of 
distinguishable photons. Assuming the incoming photons in mode 1 can be distin- 
guished from those in mode 2 (for instance, by polarization), then the two-photon 
wavefunction from Eq. [42] is 

IV'iv)la,26=fna^ t S^ t ]|0) ! (46) 



with [a^,6 3 ] = V The four possible (distinguishable) cases are: detecting an 
photon in mode 3 and a w photon in mode 4; detecting a w photon in mode 
3 and an a) photon in mode 4; detecting an a* photon in both mode 3 and mode 
4; and detecting a w photon in both mode 3 and mode 4. Thus, the second-order 
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correlation function is 



Pa,b,j(t,td) — la,2b (ii N \Ez a (t)E-(t + t d )E+(t + t d )E+ a (t) \^ N ) 



2b(^N\E- b (t)E^ a (t + t d )E+ a {t + t d )E+{t)\^ N ) 



/la, 26 



+ la,26 (^ N \E^ a (t)E^ a (t + t d )E+ a (t + t d )E+ a {t)\i> N ) la, 26 
+ la,26 (iP N \E- b (t)E-(t + t d )E+(t + t d )E+(t)\1>N) la, 26 



2la,26 



(i> N \Ez a (t)E; b (t + t d )E+(t + t d )E+{t)\i> N ) 



la. 26 



+2 la, 26 ^ N \E^{t)E^ a {t + t d )E+ a {t + t d )E+(t)\ij> N ) la, 26, (47) 

where the last line follows from symmetry. In addition, since [aj,Sj] = V i,j, the 
second term in Eq. 47 is equal to (two times) the second-order correlation function 
for a train of single-photons (Eq. 36 1 . On the other hand, the first term is 



2 la , 2b (^N\Es a (t)E- b (t + t d )E+(t + t d )E+(t)\ip N ) la , 2b 
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Writing the operators for modes 3 and 4 in terms of 1 and 2 yields 

2 la , 26 <^ V |£ 3a (i)£ 4fe (* + t d )i? 4 + b (t + id)£ 3 + aWI<Mla,26 
JV AT 
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(48) 
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(f[ a^A ^H^b^a^ (f[ S^>% m A 

\m=0 J \m=0 ) 



N N 



= 2^2 E E e -(.t-kt P )/r e -(t+t d -nt p )/r e ^ _ fc<p ) Q ( t + ^ _ ^ (49) 



Plugging this into Eq. 47 and integrating over i (again, integrals identical to Eq. 37 1, 
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Fig. 9. Theoretical two-photon second-order correlation function, for identical (blue, solid line; 
Eq. |45| and distinguishable (red, dotted line; Eq. \50\ photons. Plotted is the relative probability 
of detecting a photon in both mode 3 and mode 4 as a function of the time delay (t^) between 
the two photons for N = 20. The repetition rate shown is 20 times the natural lifetime of the 
atomic excited state (r). The decrease in the peak amplitude of the probability of joint detection 
for increasing delay time results from having a finite train of impinging photons. Interference 
between identical photons impinging on the beamsplitter results in no events at zero time delay 
= (blue, solid line). In contrast, distinguishable photons do not interfere, and thus coincident 
detections are allowed, as evidenced by the peak at td = (red, dotted line). 



we obtain the final form for the second-order correlation function for distinguishable 
photons: 



N N 



)l/r 



(50) 



k=0 n=0 



This result for distinguishable photons impinging on the beamsplitter is plotted 
against the case of identical photons (Eq. 45 1 in Fig. [9] 



4.2. Experimental interference of photons at a beamsplitter 

The interference of single photons at a beamsplitter was accomplished by using 
two trapped 174 Yb + atoms as single photon emittersP3 A single ytterbium ion is 
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confined in each of two nearly identical four-rod rf traps that are contained in in- 
dependent vacuum chambers. As described in Sec. [3j the ions are Doppler-cooled 
by cw laser light at 369.5 nm, and the presence of a single atom in each trap is con- 
firmed by imaging this fluorescence on a camera. In addition to the objective lens 
used to observe the ion fluorescence on the camera, we position a second objective 
with effective numerical aperture ~0.3 on the other side of the vacuum chamber, as 
illustrated in Fig. [10] Photons collected by this objective are coupled into a single- 
mode optical fiber and directed to a 50:50 nonpolarizing beamsplitter Q The output 
of the fiber is incident on the beamsplitter at a small angle (about 10°) to obtain 
polarization independence. The photons are detected by two PMTs, positioned at 
each exit port of the beamsplitter. Polarizers on motorized flip-mounts can be added 
to the beam path between the beamsplitter and each PMT to enable detection of 
only identically (parallel) polarized photons. Conversely, distinguishable (perpen- 
dicularly polarized) photons are detected by the addition of a A/2-waveplate to one 
of these beam paths. The arrival times of the photons are recorded by a time-to- 
digital converter (TDC) with 4 ps resolution F] 




Fig. 10. Experimental setup for the two-photon interference experiment. Fluorescence of an 
Yb+atom during Doppler cooling can be monitored by either a photon-counting photomultiplicr 
(PMT) or a camera (the grayed mirror is removable by a motorized flip-mount). The second-order 
correlation function (joint detection probability) is measured by coupling spontaneously emitted 
photons from the ion into a single-mode optical fiber, the output of which is directed to a 50:50 
nonpolarizing beamsplitter (BS). Polarizers (PBS) and a A/2-waveplate (A/2) are removable for 
the measurement of unpolarized, parallel polarized, or perpendicularly polarized photons. 

Ultrafast excitation of an Yb + atom from 2 Si/ 2 to 2 P\/2 results in the generation 
of a single spontaneously emitted photon. We use pulses of about 1 ps duration with 
central wavelength at 369.5 nm to drive the 2 S , 1 / 2 2 P\/2 transition in the atom. 

'The optical fiber is StockerYale NUV-320-K1. The fiber is single-mode at 370 nm, with a specified 

attenuation of approximately 0.1 dB/meter. 

iThe TDC is a PicoHarp 300, made by PicoQuant. 
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Since the duration of the pulse is much shorter than the excited state lifetime (8.12 
ns)p^the probability that the atom scatters two photons during a single pulse is 
approximately given by the probability of one spontaneous decay over the duration 
of the pulse: 1 - e"^ P S )A 8 - 12 ns ) « 1(T 4 . 

The 1 ps ultrafast pulses are produced by an actively mode-locked Ti: sapphire 
laser (Spectra-Physics Tsunami) operating at a central wavelength of 739 nm and 
a repetition rate of about 81 MHz. This pulsed laser is pumped by approximately 6 
W of 532 nm light (Spectra-Physics Millennia) , producing an average output power 
of 1 W at 739 nm. We measured the width of the pulses to be 1 ps by using a sim- 
ple autocorrclator, consisting of a Michelson interferometer and a lithium triborate 
(LBO) crystal with a photodetector. The pulses are passed through an electro-optic 
pulse picker that has an average extinction ratio better than 100:1 in the infrared. 
With the pulse picker synced to the intracavity acousto-optic modulator (AOM) 
of the pulse laser, we reduce the effective pulse repetition rate to about 8.1 MHz 
by allowing only one in every ten pulses to pass. A critically phase-matched LBO 
crystal is used to frequency double each pulse, resulting in approximately 80 mW 
average power at 369.5 nm. This second-harmonic light is separated from the fun- 
damental infrared by a prism. In addition to creating the wavelength needed to 
resonantly excite the ion, frequency-doubling increases the effective average extinc- 
tion ratio of the pulse picker to about 10 4 : 10 Even though the second-harmonic 
generation efficiency is only about 8%, the resulting pulses at 369.5 nm have more 
than enough power to excite an atom with unit probability. This is demonstrated in 
Fig. [II] where the scattering rate of a single atom is measured as a function of the 
incident pulse power, indicating each pulse has enough power to drive about a 3.57T 
rotation between the 2 Si/2 and 2 P\/2 levels. As the Rabi frequency is proportional 
to the square root of the incident pulse power, this means that a single pulse could 
be split into about 12 pulses, each with enough power to drive a it rotation.^ In 
the experiments presented here, each pulse is attenuated to achieve approximately 
unit excitation probability (l7r pulse). 

Before looking at the interference of two photons, we first demonstrate that the 
ultrafast excitation of a single Yb + atom is an excellent source of single photons. For 
this measurement, only the photons scattered by one of the two ions are incident on 
the beamsplitter (the other ion is blocked) . The experiment consists of a repetitive 
sequence of 10 us of cooling followed by 40 /is of excitation/measurement. During 
the cooling interval, the pulse picker is switched off so that the ion is illuminated 
solely by the cw lasers used for Doppler cooling. For the excitation/measurement 
interval the 369.5 nm cw light is switched off by an AOM, and the pulse picker is 
gated open to allow a train of 1 ps pulses, separated in time by about 124 ns, to 
sequentially excite the atom from 2 Si/ 2 to 2 P\/2- After each excitation, the atom 

k The immediately trailing pulse is not extinquished as well as the other pulses (only about 10 3 :1), 
due to the finite switching of the bias voltage applied to the pulse picker. Averaging over all pulses 
though yields an extinction ratio of approximately 10 :1. 
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Fig. 11. Resonant picosecond Rabi flopping. The scattering rate of a single 174 Yb + atom is mon- 
itored as a train of 1 ps pulses of variable power excite the atom. The Rabi frequency for driving 
the atom from 2 S]/2 to 2 P\/2 is proportional to the square-root of the incident pulse power. At 
maximum pulse energy, we achieve about a 3.57T rotation. The green curve is a fit to the function 
Asm 2 (By r PpL/2), where Ppi, is the average power of the pulsed laser incident on the ion, and 
A and B are the fit parameters. 

will spontaneously decay while emitting a single photon. Emitted photons at 369.5 
nm are collected by the objective lens, coupled into the single-mode optical fiber, 
directed onto the beamsplitter, and detected by the PMTs. The arrivals times of 
the detected photons allow us to determine the joint detection probability, shown in 
Fig. [12] The functional form of the experimental joint detection probability matches 
the theorectical single-photon source illustrated previously (Fig. [8]), with the lack 
of a peak at delay time td = indicating that at most a single photon is emitted 
after each ultrafast excitation pulsed Dark counts on the PMTs result in a small 
background contribution at all delay times; coupling the light from the atom into a 
single- mode fiber highly suppresses the contribution of background scattered light. 

We observe the quantum two-photon interference effect by simultaneously ex- 
citing both trapped Yb + atoms, and combining the emitted photons on the beam- 
splitterQ The joint probability of detection for identical photons is measured by 
placing both polarizers in the beam path, while detection of distinguishable pho- 
tons is evaluated by including the A/2-waveplate in one of the two paths. As was 
seen theoretically in Fig. [9j for identical photons we expect no detections at time 

'Of course, simultaneously excitation is not essential; the emitted photons just need to arrive at 
the beamsplitter simultaneously. As such, the path length of the excitation pulses are adjusted so 
that emitted photons from each ion arrive at the beamsplitter within 100 ps of each other. Since 
the path lengths between each ion and the beamsplitter are approximately equal, this means the 
ions are also excited at about the same time. 
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Fig. 12. Experimental single-photon second-order correlation function. The data is shown with 1 
ns binning, and was integrated over 3 hours. The lack of a peak at time delay = is indicative 
of a source of single photons, confirming that at most one photon is emitted by the atom following 
ultrafast excitationEZl 

delay ta = 0, while for distinguishable photons we expect a probability of joint 
detection half as large as that from adjacent excitations. The data shown in Fig. [T3| 
demonstrates this quantum two-photon interference effect. In the measurement of 
parallel polarized photons, the residual counts at time delay t& — result from 
both dark counts on the PMTs, and imperfect spatial mode overlap of the photons 
on the beamsplitter. The data shown in Fig. [13] corresponds to an interferometer 
contrast of approximately 95%. 87 [^]The spatial filtering afforded by coupling the 
spontaneously emitted photons into single-mode fibers was essential to achieving 
this two-photon interference. 

5. Remote Quantum Gate 

The standard model of quantum computation is the quantum circuit model, where 
a series of gate operations are performed and at the output the qubits are measured 
as a final step. Deterministic, two-qubit gates that utilize the common modes of 
motion of trapped atomic ions have already been achieved with 99.3% fidelity.^ 
Effort in this area is focused on scaling the system to larger numbers of qubits, as 
can be done, for instance, by constructing complex trap arrays where ions can be 
stored and operated on at various intervals and locations.^ Crucial to the success 
of this methodology will be understanding the motional decoherence that plagues 
all ion trapsEEZElEll 

Another approach to scalable quantum computation is the cluster-state (or one- 
way) quantum computing model, where the measurement process itself is an integral 
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m More recently, the contrast of our interferometer has reached > 98%, as measured with laser 
light (Sec. 16.3.1 1. 
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Fig. 13. Experimental two-photon second-order correlation function for: (a) distinguishable (per- 
pendicularly polarized) photons; and (b) identical (parallel polarized) photons.EZlBoth sets of data 
are shown with 1 ns binning, and each data set was integrated over about 4 hours. As expected 
from the theorectical calculations of Sec. |4.1| two identical photons incident on the beamsplitter 
always exit by the same port, resulting in suppression of joint detections at time delay td = in 
(b). 



part of the gate operation.^ In this scheme, a large entangled state is prepared first. 
The computation is then executed by successive measurements of the qubits, with 
classical feed-forward of prior measurement results determining the single-qubit 
gates to be performed on subsequent sets of qubits before their measurement. This 
approach has been proven to be equivalent to the circuit modelj^ and allows for 
a new class of quantum gates to be considered for scalable quantum computation: 
non-unitary, measurement-based gates P 

Basic quantum operations applicable to the cluster-state model of quantum 
computation have been carried out using entangled photons However, the pho- 
tonic cluster states in these experiments were generated by spontaneous parametric 
down-conversion, relied on post-selection, and did not utilize a quantum memory, 
making these implementations difficult to scale to larger systems.^ 

We review the implementation of a probabilistic, heralded quantum gate be- 
tween two remote quantum memories. A full set of input states is used to evaluate 
the operation of the gate, resulting in an average fidelity of 89(2)%PEven though 
the gate is intrinsically probabilistic, the incorporation of a quantum memory allows 
this system to be efficiently scaled, potentially allowing for the generation of the 
large cluster-states necessary for one-way quantum computationJH2ME3 While in 
this demonstration the success probability is only about 2.2 x 1CP 8 , it should be pos- 
sible to improve the success rate for practical implementation. The photon-mediated 
approach to scalable quantum computing with trapped atomic ions has the advan- 
tage that complex trap arrays may not be required, the operation is insensitive to 
motional decoherence, and the gate can be operated on qubits separated by a large 
distance. Alternatively, spatially separated quantum registers of Coulomb-coupled 
ions could be linked using this heralded gate to establish a quantum network. More- 
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over, since the operation is mediated by photons, this gate could be applied to any 
system of optically active qubits, such as neutral atoms, ions, nitrogen- vacancy cen- 
ters, or quantum dots. Indeed, hybrid systems are also envisioned, where disparate 
quantum systems are connected via this photon-mediated process to exploit the 
advantages of each individual quantum system. 



5.1. Ion-photon entanglement 

The implementation of the gate begins by confining and cooling two 171 Yb + atoms 
in two vacuum chambers, separated by a distance of about one meter. The exper- 



imental setup, shown in Fig. 14 is similar to that employed for the two-photon 



interference experiments of the previous section. In this case, though, the 171 iso- 
tope of Yb + is used, and the ions are subjected to an external magnetic field of 
about 5.2 gauss that is aligned perpendicular to both the observation axis and the 
impinging light from the picosecond pulsed laser. Also, the polarizers arc present 
in the setup, whereas the optional A/2-waveplate used previously is removed from 
the beam path. 




Fig. 14. Experimental setup for the heralded quantum gate. Spontaneously emitted 7r-polarized 
photons arc coupled in a single-mode fiber and directed to interfere on a 50:50 non-polarizing beam- 
splitter (BS). Coincident detection of two photons by photon-counting photomultipliers (PMTs) 
announces the success of the gate between the two ions. Polarizers (PBSs) are used to filter the 
photons so that only 7r-polarized photons are detected. The state of each ion is measured by 
state-dependent fluorescence, detected by a PMT on the opposite side of the vacuum chamber. 



After Doppler cooling, each ion is initialized by optically pumping to the state 
|0) using a 1 /is pulse of 369.5 nm light. Each of the two ions, ion a and ion b, 
are then prepared in a particular quantum state by resonant microwave radiation 
at 12.6 GHz. The microwaves are applied directly to one of the trap electrodes in 
each vacuum chamber with controlled phase and duration (0-16 /is). Thus, after 



Quantum Logic Between Distant Trapped Ions 31 



application of the microwave radiation we write the state of the atoms as 

\i>) a = a\0) a + p\l) a 

|^>6 = 7|0> 6 + <5|l> b (51) 

with \a\ 2 + \{3\ 2 = 1 and | — ^ | 2 + \5\ 2 = 1. For clarity, we will assume ideal state 
evolution throughout this discussion. 

After state preparation, each ion is excited with near-unit probability to the 
2 P\/2 level by an ultrafast (1 ps) laser pulse having linear polarization aligned 
parallel to the magnetic field (7r-polarized) and central wavelength at 369.5 nm. 
Due to the polarization of the pulse and atomic selection rules, the broadband 
pulse coherently drives |0) to |0') := 2 P 1/2 \F = 1, m F = 0) and |1) to |1') := 



2 P\/2\F — 0,mp = 0), as illustrated in Fig. 15 'a). Since the bandwidth of the 1 



ps pulse is approximately 300 GHz, we are able to drive both of these transitions 
simultaneously.^1 The 100 THz fine structure splitting of the 2 P levels ensures that 
the coupling to 2 P^,/2 is negligible. 
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Fig. 15. Pulsed excitation with 7r-polarized light to generate ion-photon entanglement, (a) Due 
to the atomic selection rules, 7r-polarized light drives |0) to |0'} and |1) to |1') only. The pulse 
bandwidth of about 300 GHz allows both of these transitions to be driven simultaneously with 
near-unit excitation probability, (b) After excitation, the ion spontaneously decays back to 2 S\i2 
and emits a single photon at 369.5 nm. If we consider only 7r-polarized photons, then the fre- 
quency of the emitted photon is entangled with the internal electronic state of the atom, with the 
separation of the different frequency modes equal to the sum of the 2 Si/2 an d 2 ^ > i/2 hyperfine 
splittings (yblue. ~ "red = 14.7 GHz). Polarizers (PBS) are used to filter out the cr-polarized light. 



As seen in Sec. |4j each ion spontaneously emits only a single photon at 369.5 
nm while returning to the 2 <S I 1 / 2 ground state. The emitted photons at 369.5 nm 
can each be collected along a direction perpendicular to the quantization axis by an 
objective lens with effective numerical aperture ~0.3 and coupled into a single-mode 
fiber. Observation along this direction allows for polarization filtering of the emitted 
photons because those produced by 7r and a transitions appear as orthogonally 
polarized.^ Considering only 7r decays results in each ion being entangled with the 
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frequency of its emitted photon such that 



\lp)a = a\0) a \l/ b lue)a + P\ l)a Wred) a 
\lp)b = l\fybWblue)b + S\l)b\v red ) b , 



(52) 



where \vuue) and \v T ed) are the single photon states having well-resolved frequencies 
Vbiue and v rec i, respectively. The bandwidth of each state is determined by the nat- 
ural linewidth of the excited state to be 1/(2itt) = 19.6 MHz.^The two frequency 
states are separated by the sum of the 2 Si/2 and 2 P\/2 hyperfine splittings, so that 
Vbiue — v re d = 14.7 GHz. The output of the single-mode fibers is then directed to 
interfere at a 50:50 nonpolarizing beamsplitter. 



5.2. Ion-ion entanglement 

Although the action of the beamsplitter on identical and distinguishable photons 
was reviewed in Sec. |4j the case of superposition states requires additional treat- 
ment. In order to evaluate the action of the beamsplitter in the present scenario, 
the state of the two ion-photon entangled systems from the previous section is ex- 
pressed in the Bell state basis of the photons. The total state of the two ion-photon 



systems is simply the product of the states in Eq. 52 

l*> = IV>}a® \i/>)b 

= 07|0)«|0) 6 | \Vblue)b + a<5|0) a |l) 6 | ^blue) a\^red) b 

+f3l\l)a\0)bWred)aWblue)b + @6\ l) a | 1) bWred) a Wred) b 

= a 7 |0) a |0) h -^ {\ct> + ) ph + \<f>-) P h) + aS\0) a \l) b ^= (\1>+) ph + \rP~) ph ) 
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+ \-4> + ) P h M0).|l>6 + /3 7 |1) Q |0) 6 ) + \i3-)ph (a5|0)„|l) 6 - fr\l) a \0) b )} (53) 

The action of the beamsplitter can then be determined for the four possible Bell 
states. Assuming photon a enters the beamsplitter via port 1 and photon b enters 
the beamsplitter by port 2, and associating the creation operators a) and w with 
frequency modes \v re d) and \vuue)i respectively, we then find 
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for the \(j> ± ) p h state, where in the last line the notation shows the number of photons 
of a particular frequency in the spatial mode 3 or 4. Similarly, for the \ip + ) p h and 
\'4 ) ~}ph photon states, we find 
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The calculations above show the striking result that only the antisymmetric Bell 
state, \ifj~) p i u produces a single photon in both exit ports of the beamsplitterP3 
This crucial point is the final step in the quantum gate. 

We use the joint detection of two photons at the exit ports of the beamsplitter 
to herald the success of the quantum gate. By the above derivation, a coincident 



detection signals that the photons were measured in the \ip ) p h state. From Eq. 53 
the coincident detection thereby projects the ions into the state 



(a5|0) Q |l> 6 -/? 7 |l) a |0) b ), 



(57) 



where the factor i9 = (|a| 2 |<5| 2 + |/3| 2 |7| 2 )/2 in front is the usual renormalization 
term present after a measurement.^ In the case |a| = |/3| — \y\ = \S\ = l/y/2, 
the ions arc left in a maximally entangled state. In terms of operators, this gate is 
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written as 

^ 3 °(<7 V 6 - o$4), (58) 

where <r] is the jth Pauli operator acting on the ith qubitP^ 

In contrast to our earlier demonstrations of remote ion entanglement the 
initial state amplitudes are preserved by the heralded quantum gate and determine 
the form of the final ion-ion entangled state. This is a defining feature of a gate 
operation, and is essential to establishing entanglement between more than two 
qubits. However, the quantum gate presented here is not unitary. Indeed, for certain 
initial states (e.g. a = 7 = 1), a coincident detection should never occur. While this 
behavior voids its application to the quantum circuit model, cluster states can still 
be generated by having all qubits initially in a superposition state. In this case, the 
gate succeeds with nonvanishing probability, and scales favorably.^ 



Table 2. Experimental characterization of the heralded quantum gate. Listed are the input states of the two ions, the 
expected (ideal) output state after operation of the gate, the measurement performed, the number of heralding events, the 
obtained fidelity, and the measured and ideal probability for two photons to be in the antisymmetric Bell state. The fidelity 
of the output states is obtained by parity measurements in the appropriate bases. Here we have defined the parity P xy as the 
difference of the probabilities to find the two ions in the same state and in opposite states when ion a, b is measured in the x, y 
basis, respectively. The other parity values are defined similarly (/ indicates independence to that qubit measurement). From 
these results we calculate an average fidelity T = 0.89(2)P1 The success probability of the gate is P ga te = J? x (8.5) X 10 8 . 



measured theory 



input state expected output measurement events fidelity 



(|0>a + |l>a)Cgl(|0>6 + |l>6) |0>a|l>6 - |l>a|0> 6 | (1 - P xx - Pyy - Pzz ) 210 0.89(2) 0.26(1) 1/1 

(|0)a+i|l)o)® (|0) b +|l)i,) |0)a|l>6-i|l)a|0) i , | (1 - P x y + Pyx - Pzz ) 179 0.86(2) 0.26(1) 1/4 

(|0) a -|l)a)® (|0)(, + |l)f,) |0) a |l)6 + |l) a |0)i, | (1 + P xx + P yy - P zz ) 178 0.85(1) 0.22(2) 1/4 

(|0) a -i|l) a )®(|0)(, + |l)i,) |0)a|l> 6 + i|l>a|0> b 1(1 + Pxy - Pyx - Pzz) 188 0.81(2) 0.27(2) 1/4 

([Q)o + |l)a)®[l)6 |0)a|l)i, -A1 + Plz- Pzl ~ Pzz) 42 0.86(5) 0.24(4) 1/4 

|0>a® (|0> 6 + |1>6) |0)a|l)f, | (1 + P Iz - P zI - P M ) 52 0.90(4) 0.20(3) 1/4 

|0> a (g)|l>6 |0>a|l>6 \ ( 1 + P Jz - P z i - Pzz ) 48 0.98(2) 0.39(6) 1/2 

|0} a ®|0} 6 65 0.04(1) 



5.3. Gate evaluation 

We verify the operation of the gate by implementing it on a variety of input states. 
The gate is characterized by determining the fidelity of the resulting output state 
with respect to the ideal case described by Eq. [57] The fidelity is defined as the 
overlap of the measured state with the ideal state 

? '= {i>ideal\PWideal) , (59) 

where we have written our measured density matrix as p, and the ideal output 
state as \ipideai}- The various input states used to characterize the quantum gate, 
together with the measurements made and the resulting fidelities for each, are 
given in Table [2] As measurement of each ion is accomplished using the afore- 
mentioned state fluorescence technique, measurement in the remaining two bases 



Quantum Logic Between Distant Trapped Ions 35 



requires an additional microwave pulse before detection. We define the rotation 
{Ry(w/2), R x (tt/2),R(0)} before detection to correspond to measurement in the 
basis {x,y, z}. Overall, we attain an average fidelity of 89(2)%. 94 

To evaluate the gate, we do not make measurements of all combinations of pos- 
sible input states around the equator of the Bloch sphere, as input states that differ 
only by a global microwave phase rotation on the equator of the Bloch sphere are 
indistinguishable in the lab. It is only the phase difference between the states that 
matters experimentally, where the microwave oscillator is used as a reference. As an 
example, consider the two possible combinations of input states l/v2(|0) a + |l) a ) 
and 1/V2(|0} 6 + i\l)b) versus l/y/2(\0) a + and l/v^(-|0) 6 + |1) 6 ). Experi- 

mentally, the preparation of both of these pairs of states is the same: application 
of a microwave 7r/2-pulse to ion a, and application of a microwave 7r/2-pulse to ion 
b that is +7r/2 out of phase with the pulse that was applied to ion a. 

In addition to the fidelity measurements of Table [2] we also assess the generation 
of the maximally entangled antisymmetric Bell state \ipideai) — l/v / 2(|0) a |l)b — 



|l) a |0)b) by full state tomography. The resulting density matrix shown in Fig. 16 



is obtained using a maximum likelihood method P^Sl From this density matrix we 
calculate an entangled state fidelity of T — 0.87(2), a concurrence of C = 0.77(4) 
and an entanglement of formation Ep — 0.69(6)1^1 




Fig. 16. State tomography of p for the case \ipideal > = (|0>a|l>6 - |l>a|0> 6 )/%/2. The figure shows 
the real (a) and imaginary (b) elements of the reconstructed density matrix. The density matrix 
was obtained with a maximum likelihood method from 601 events measured in 9 different bases. 
From this density matrix we calculate an entangled state fidelity of T = 0.87(2), a concurrence of 
C = 0.77(4) and an entanglement of formation Ep = 0.69(6)S 

The observed fidelity is consistent with known experimental errors. The pri- 
mary contributions to the error are imperfect state detection (3%), spatial mode 
mismatch on the beamsplitter (6%), and detection of cr-polarized light due to 
the finite solid angle of collection and misalignment of the magnetic field (<2%). 
Other sources, including imperfect state preparation, pulsed excitation to the wrong 
atomic state, dark counts of the PMT leading to false coincidence events, and mul- 
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tiple excitation due to pulsed laser light leakage, are each estimated to contribute 
to the overall error by much less than 1%. Micromotion at the rf-drive frequency 
of the ion trap, which alters the spectrum of the emitted photons and can degrade 
the quantum interference, is expected to contribute to the overall error by less than 
1%. 

5.4. Success probability 

The quantum gate is a heralded, probabilistic process. The net probability for 
coincident detection of two emitted photons is given by 

Pgate = \p^T fiher T optics ^ Afl/Air)} 2 « x (8.5 x ltr 8 ) , (60) 

where p n — 0.5 is the fraction of photons with the correct polarization (half are 
filtered out as being produced by a decays); 77 = 0.15 is the quantum efficiency of 
each PMT; T/j(, er = 0.2 is the coupling into, and transmission through, the single- 
mode optical fiber; T op ti CS = 0.95 is the transmission through the other optical 
components; £ = 1—0.005 = 0.995, where 0.005 is the branching ratio into the 2 D 3 / 2 
level; Ail/4ir — 0.02 is the solid angle of light collection; and •§ is the normalization 
factor from Eq. [57] whose value depends on the initial amplitudes of the ion qubit 
states. 

In the current setup, the attempt rate was limited to about 75 kHz, due to the 
time required for the state preparation microwave pulse. This resulted in about one 
successful gate operation every 12 minutes. However, the expression for P ga te reveals 
multiple ways to substantially increase the success rate. The most dramatic increase 
would be achieved by increasing the effective solid angle of collection, which, for 
instance, could be accomplished by surrounding each ion with an optical cavity. 
Although improvements that increase the success probability of the gate operation 
can enhance scalability, even with a low success probability this gate is still scalable 
to more complex systems.^ 

6. Quantum Teleportation 

Quantum teleportation is the faithful transfer of quantum states between systems 
that relics on the prior establishment of entanglement, but uses only classical com- 
munication during the transmission.^ Teleportation, where quantum information 
is transferred between two disparate locations without traversing the space be- 
tween the systems, is a stark realization of the counter-intuitive aspects of quan- 
tum physics. In addition, the ability to tclcport quantum information is an essen- 
tial ingredient for the long-distance quantum communication afforded by quantum 
repeaters^ and may be a vital component to quantum computation.^^ 

Previous demonstrations of quantum teleportation have been accomplished sep- 
arately using optical systems over long distances) 102 ! 103 ! 104 ! 18 ! 105 ! and quantum 
memories in close proximity) 106 ! 107 !!^ In this section, we review the implemen- 
tation of a heralded teleportation protocol where the advantages from both optical 



Quantum Logic Between Distant Trapped Ions 37 



systems and quantum memories are combined to teleport quantum states between 
two trapped ytterbium ion qubits over a distance of about one meter. The protocol 
is evaluated on a complete set of mutually unbiased basis states, resulting in an 
average teleportation fidelity of 90(2)%P^The execution of the heralded quantum 
gate in this basic quantum algorithm elucidates the applicability of this scheme for 
quantum networks. 



6.1. Teleportation protocol 

The experimental setup for the teleportation protocol is identical to the one used 



to demonstrate the quantum gate (Fig. 14). A single 171 Yb + atom is confined and 



cooled in each of two nearly-identical traps separated by a distance of about one 
meter. The atoms are initialized in the state |0) by optical pumping via a 1 (is pulse 
of 369.5 nm light. We then prepare the ions a and b in the states 

\f) a = a\Q) a +l3\l) a (61) 

|V)b = ^(|0) b + |l) fe ) 

by applying a resonant microwave pulse of controlled phase and duration (0-16 /is) 
directly to one of the trap electrodes, where \a\ 2 + \(3\ 2 = 1. The quantum state 
written into ion a is the information we seek to teleport. While in the present case 
the amplitudes are determined by the applied microwave radiation, in principle a 
and (3 could be unknown. 

Following state preparation, each atom is excited by a 1 ps pulse of 7r-polarized 
light at 369.5 nm, and this broadband excitation coherently drives the population 
in the hypcrfine levels of the 2 S\/2 state to complementary levels in the 2 P\/2 state, 



as shown previously (Fig. 15h . As each atom decays back to the 2 S\/ 2 levels, it emits 
a single photon at 369.5 nm. Considering only ir decays results in the frequency of 
the emitted photon becoming entangled with the electronic state of the ion, such 
that: 

\lp)a = Qi|0) a |z/b/ ue ) a + f3\l) aWred) a 

\ij)) b = -j= (|0) 6 |^„ e ) 6 + \l) b \u red ) b ) . (62) 

Spontaneously emitted photons are collected by an objective lens, coupled into a 
single-mode fiber, and directed to interfere at a 50:50 non-polarizing beamsplitter. 
As shown in Sec. [5] due to the quantum interference of the photons at the beam- 
splitter, a simultaneous detection at both output ports of the beamsplitter occurs 
only if the photons are in the \i/j~) p h state, so that the action of this heralded 
quantum gate is to project the ions into the entangled state: 

l^hons = ph{i>~\ (Wa ® \lp)b) = a\0) a \l) b - /3|l) a |0)fc. (63) 

Even though the gate is a probabilistic process, we do not require postselection 
because the coincident detection of the two photons is a heralding event that an- 
nounces the success of the quantum gate. In the teleportation protocol d = 1/4 in 
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Eq. 60 so that the overall success probability for any a and is about 2.2 x 10 , 
limited by the efficiency of collecting and detecting both spontaneously emitted 
photons. Therefore, the previous steps (state preparation and pulsed excitation) 
are repeated at a rate of 40 to 75 kHz, including intermittent cooling, until the gate 
operation is successful (once every 12 minutes, on average). 

After the success of the quantum gate has been confirmed by the heralding 
event, ion a is subjected to another pulse of microwaves to execute the rotation 
Ry(n/2). The microwave rotation transforms the state given in Eq. 63 to: 



\i>) lons = ^ (|0) a + |l) a ) |1) 6 - -L (_| ) o + |l) a ) |0) b 

= i (a\l) b + 0\O) b ) \0) a + i (a\l) b - 0\O) b ) \l) a . (64) 

We then measure ion a using the state dependent fluorescence technique dis- 
cussed earlier. As is apparent from Eq. |64| measuring ion a projects ion b into one 
of the two states: 

If measured |0) a => \tj)) b = a\l) b + 0\O) b 

If measured |1) =*> \ip) b = a\l) b - 0\O) b . (65) 

The result of the measurement on ion a is relayed through a classical com- 
munication channel and used to determine the necessary phase of a conditional 
microwave 7r pulse applied to ion b to recover the state initially written to ion a. If 
\0) a is measured, the rotation R x (n) is applied to ion b so that: 

Rx(n) (a\l) b + 0\O) b )=a\O) b +0\l) b . 

On the other hand, if |l) a is measured, the rotation Ry(jr) is applied to ion 6, 
yielding 

R y (n) (a|l) 6 -/3|0) 6 ) = a|0) h + /3|l) i) , 

where we have ignored any global phase factors. Thus, the state amplitudes initially 
written to ion a are now stored in ion b, which completes the teleportation of the 
quantum state between the two distant matter qubits. 



6.2. Experimental evaluation 

We execute the teleportation protocol on a set of six mutually unbiased basis states 

I^W e {iA/2(|o>+|i», iA/2(|o>-|i», i/V2(|o>+i|i», iA/2(|o)-i|i», |o), |i)}, 

and evaluate the process by performing state tomography on each teleported state. 
A single-qubit density matrix can be reconstructed by measuring the state in 
three mutually unbiased measurement bases. Of course, measurement of the ion 
occurs via the aforementioned state fluorescence technique, and therefore only 
distinguishes between |0) and |1) (z-basis); two states such as l/y/2(\0) + |1)) 
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and 1 / v^2(|0) — |1)) (x-basis) are not distinguishable by fluorescence alone. Mea- 
surement in the remaining two bases requires an additional microwave pulse be- 
fore detection. As before, rotation {R y (n/2), R x (ir/2), R(0)} before detection to 
correspond to measurement in the basis {x,?/, z}. These measurements allow re- 
construction of the single-qubit density matrix, p, for each teleported state us- 
ing a simple analytical expressionJSSl with the results shown in Fig. 17 The fi- 
delity of the teleportation protocol, defined as the overlap of the ideal and mea- 
sured density matrices T = tr(p idea ip) = idea i (i/)\p\ip) ideal, for this set of states 
is measured to be T = {0.91(3), 0.88(4), 0.92(4), 0.91(4), 0.93(4), 0.88(4)}. This 
yields an average teleportation fidelity of T = 0.90(2) J^H The reconstructed den- 
sity matrices also facilitate characterization of the protocol by quantum process 
tomography.^ Using a maximum likelihood method)*^ we determine the process 
fidelity T process = tr(\ideaix) = 0.84(2), which is consistent with the aforemen- 
tioned average fidelity of the teleportation protocol!^] 
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Fig. 17. Tomography of the teleported quantum states. The reconstructed density matrices for 
the six unbiased basis states teleported from ion a to ion b: (a) \ipideal) = 1 / \/2( |0) + 1 1} ) teleported 
with fidelity T = 0.91(3), (b) \ip idea i) = l/%/2(|0> - |1» teleported with fidelity T = 0.88(4), (c) 
\f ideal) = Vv^(|0) + i|l>) teleported with fidelity T = 0.92(4), (d) \ij> ideal ) = l/%/2(|0> - i|l)) 
teleported with fidelity T = 0.91(4), (e) \ip ldea i) = |0) teleported with fidelity T = 0.93(4), and 
(f) \1pideal) = |1) teleported with fidelity T = 0.88(4). These measurements yield an average 
teleportation fidelity T = 0.90(2). The data shown is from 1285 events recorded over 253 hours. 
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6.3. Error analysis 

The experimental errors in the teleportation protocol differ slightly from the im- 
plementation of the quantum gate. The primary sources of error that reduce the 
average fidelity are photon mode mismatch at the 50:50 beamsplitter (4%), im- 
perfect state detection (3.5%), and polarization mixing resulting from the nonzero 
numerical aperture of the objective lens and from misalignment with respect to the 
magnetic field (2%). Other sources, including incomplete state preparation, pulsed 
excitation to the wrong atomic state, dark counts of the PMT leading to false coin- 
cidence events, photon polarization rotation while traversing the optical fiber, and 
multiple excitation resulting from pulsed laser light leakage, are each expected to 
contribute to the error by much less than 1%. Residual micromotion at the rf-drive 
frequency of the ion trap, which alters the spectrum of the emitted photons and 
degrades the quantum interference, reduces the average fidelity by less than 1%. 
Detailed calculations of the two most prominent error contributions are presented 
below. 



6.3.1. Photon spatial mode-mismatch at the beamsplitter 

We determine the effect of imperfect spatial mode-overlap of two photons incident 
on the beamsplitter by utilizing the formalism presented in Sec. |4j In this case, 
only two photons impinging on the beamsplitter are considered, rather than the 
train of photons described previously. Given that photons from different excitation 
pulses are well-separated in time, this simplification can be made without loss of 
generality. 

To determine the decrease in fidelity resulting from spatial mode-mismatch of 
the photons at the beamsplitter, we calculate the effect of the electric field operators 
on the four possible input states \v re d)aWred)b, Wred)a\vblue)b, Wblue) a\vred)b, and 
\vbiue)a\vuue)b- As an attempt to unclutter the expressions, henceforth the states 
will not be labeled with a or b unless there are possible ambiguities. Instead, we will 
write, e.g. \v re d)aWblue)b — \RB), where the ordering of the frequencies determine 
to which subsystem they apply (e.g. in this case "R" indicates a "red" photon as 
part of system a) . 

The four possible ways to get a photon in each exit port of the beamsplitter are 
E+E+(t d ), E+E+(t d ), E+ a E+(t d ), E+ a E+ a (t d ), (66) 



32 



with w and a) as the 



where the electric field operators are defined as in Eq. 
creation operators for a \vuue) = \B) and \v red ) = \R) photon, respectively (as in 
Sec. [5]). In the subsequent calculation, we suppress the exponential decay factor, 
the Heaviside step function, and the normalization factors of Eq. [32] The present 
concern is the resulting density matrix if a coincident detection of two photons 
occurs, and in this derivation these other factors may be regarded as part of a 
normalization factor. 
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There are then 16 terms to calcuate; fortunately, many of them have an identical 
form. 



1 



t At I 



-£i&(t«i)ai&2 + 66(id)(a 2 )"J 4410) 

66(^)1 10) 



(67) 



This will also be the result for E+ b E+(t d )\BB) , as can be seen by switching a with 
b, and \R) with \B). 
Next, we find 



1 



E+E+(t d )\RR) = - Ci&i + 6&a -fi(td)&i + 6(*d)&2 



44 10) = 



0, (68) 



because ba^O) = 0. Of course, the term will vanish whenever there are more anni- 
hilation operators of a specific mode than creation operators. As such, the following 
terms also vanish: 



= E+ a Ei a {t d )\BB) = E+ a E+(t d )\BB) = E+E+(t d )\RR) 
= Et b Et a {t d )\BB) = E+E+(t d )\RR) = E+ a E+ a (t d )\RB) 
= E+ a E+ a (t d )\BR) = E+E+(t d )\BR) = E+E+(t d )\RB). 



(69) 



Of the other combinations we need to calculate, we find 

1 



E+ a E+(t d )\RB) = | [&ai+6a 2 ] \-^{t d )h + 6(f d )6 2 J a[bl\0) 



1 r 
2 



£i£i(id)ai&i - S,2^i(t d )a 2 bi 



The same answer will result from E+ b E+ a (t d )\BR). Finally, 



(70) 



E+ a E+(t d )\BR) = £ Kifii + £ 2 a 2 ] -€i(t d )6i + &(t d )&2 &J4|0) 



= - g 6fi (*<*)• 



(71) 



The result of Eq. 71 will also be the answer obtained for the term E^ b EJ a {td)\RB) 
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The terms calculated above can be used to analyze the two ion-photon systems. 
As seen previously, the system after excitation and spontaneous emission of a 7T- 
polarized photon is: 

|V>fem,p& = 4f [a\00)\BB) + P\10)\RB) + a\01)\BR) + 0\11)\RR)} ■ (72) 

To determine the density matrix of the two ions after a coincident detection of two 
photons, we use terminology similar to that in Ref. EZ1 to define the joint electric 
field measurement operator as 

M pKjk = E+E+ k (t d ), (73) 

where j, k € {a, b}. The joint detection probability from Sec. |4]is then 

Pj = p ,ME^ k (t d )E^E+E+ k (t d m ph 
= ph W M lh,jk M Ph,jk \i>)ph 

= HMl Kjk M ph j h \^) ph {^\), (74) 

where in the last two lines the definition of the electric field measurement operator 
and the usual properties of the trace were usedP^ 

In this formalism, the density matrix of the two ions will be given by a partial 
trace over the photon states of the electric field operators acting on the two-ion- 
photon system: 



tr 



ph 



Yl M lh,jk M ph,3k\lp)ion,ph( 



3,k 



\ [\a\ 2 \00){0O\(BB\Ml hijk M phJk \BB) 
+a/r|00)(10|(5i?|M^ j ^ 

+a(3* |00) (11| (BB\Ml hJk M phtjk \RR) + 0a*\lO) (00| (RB\M^ hjk M ph ^ k \BB) 
+ |/3| 2 |10)(10|(i?B|A^ i ^ 

+ |/3| 2 |10>(ll|( J RS|A/; /lijfc A^ p/l , jfc | J Ri?> + laplOlXOOK^I^^Af^^-fclSi?) 
+ap*\01)(W\(BR\M p - hJk M ph , jk \RB) + \a\ 2 \01)(01\(BR\Ml hdk M ph , jk \BR) 

+af3*\Ql){ll\(BR\Ml hdk M phJk \RR)+(3a*W 
+ \f3\ 2 \ll)(10\(RR\Ml hjk M ph , jk \RB)+Pa^ 

+ \(3\ 2 \ll)(ll\(RR\Ml hjk M ph , jk \RR)} , (75) 



where in the last several lines, the sum ■ k is implicit. Using the results for the 
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action of the different measurement operators on the states yields: 

Pious = I [|a| 2 |00)(00||ei6(M-6ei(M| 2 + + + + 

+ |/3| 2 |10)(10|(|^ 2 (^)| 2 + |^i(M| 2 ) 

+/3a*|10)<01| (-£8(*.i)6fi(*«0 -ZMtdmt(td)) +0 

+o+a/3*|oi><io| (-ei6(t d )ss(*«j) -££(*<06&(*<0) 

|a| 2 |oi)(oi| (|&fc(t d )| a + |66(*d)| 2 ) + + 

+0 + 0+ |/3| 2 |ll)(ll||ei6(^)-6a(i<l)| 2 ] ■ (76) 

Following a coincident detection in the teleportation protocol, ion a is rotated 
with a R y (7r/2) microwave pulse, and then measured. The ion-ion density matrix 
after rotation of ion a is: 

(Ry(*/2)) a ft ons (4( 7r / 2 ))a 

= ^ [M'lfc&fa) - 6a(MI 2 (|0) + |i»|o)«o| + (i|)(o| 
+I/3I 2 (I66(^)| 2 + |6£i(t d )| 2 ) (-|o> + |i»|o)(-(o| + (i|)(o| 

x(-|0) + |l»|0)((0| + (l|)(l| 
x(|0) + |l))|l)(-(0| + (l|)(0| 

+ M 2 {\tei(t d )\ 2 + |£i6(*«0l 2 ) (|0) + |1»|1)((0| + (1|)(1| 

\/3?\€&(t d ) - tei(t d )\ 2 (-\a) + |i»|i)(-(o| + (i|)(i|] . (77) 

Analogous to the electric field measurement operator, we define an operator for 
the measurement of the quantum state of the ion. In this calculation, we assume 
ideal measurement of the ion. The influence of imperfect state detection on the 
fidelity is calculated in the next sectionj^jDcfine the (ideal) measurement operators 
for the ion state as^ 

M oj = 10)^01, My = ll^ll, (78) 

where here j S {a, b} denotes operation on the jth ion. Since with the final rotation 
on ion b, conditioned on the measurement of ion a, the states will be equivalent, 
only the case of |0) a is presented below. Measurement of ion a is defined similarly 
to the case for the photons, where the density matrix of ion b is given by the partial 



n As both are independently small, it is a good approximation to consider each separately. 
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trace of a on the measurment operator and the quantum state of the ions. 
Pb = tr a (MlMoa (Ry(n/2)) a p ions (i$ (tt/2)) o ) 

= y & [{\ a \ 2 \bUU) 66 (Ml 2 + W {\ZMU)\ 2 + \tei(t d )\ 2 )) |o>(o| 
(-m(U)^i(t d ) - ZMmuKu)) |o)(i| 

-af3* (-£ih(t d )&&(t d ) - &&(U)tei(t d )) |1)(0| 

+ (M 2 (\tei(t d )\ 2 + \ZMt d )\ 2 ) + |/3| 2 |66(*c0 - 6rfi(t«0l 2 ) 1 1> <i |] . (79) 

Recall that in the teleportation protocol if |0) a is measured then ion b has a R x (tt) 
microwave pulse applied to it. Since this just flips the state of the ion, the final 
density matrix for ion b is given by: 

p" = y& km 2 {mi{u)\ 2 + \te 2 (t d )\ 2 ) + \(3\ 2 \^ 2 (t d ) tei(t d )\ 2 ) io)(oi 

-a/3* (-£ib(t d )&&{td) - |0)(1| 
-/3a* (-m(td)aUtd) - Si&t d )££(t d )) |1)(0| 

+ (M 2 |£i6(M - tei(td)\ 2 + |/3| 2 {\ZMt d )\ 2 + |^i(*d)| 2 )) |1>(1|] • (80) 

Evaluating this expression requires determination of the spatial mode factors £ 
in terms of a measured quantity, such as the visibility of the interferometer. The 
visibility is defined as 

v =j — YT~ ' (81) 

1 max \ ± mm 

where I max and I m i n are the maximum and minimum intensities of the incident 
light. If two fields are incident, with amplitudes Ei and E 2 , then I max — \Ei + E 2 \ 2 
and I m in = \E\ — ■ On the other hand, the intensity of a single field is simply 
/i = |£i| 2 . 

Suppose a function e(r) describes the amplitude of the incident light as a func- 
tion of position. However, the detector measures the intensity / incident over the 
area of the detector. Assuming the detector area completely covers the incident 
mode, then the registered intensity will be the integral over space: 

\e(r)\ 2 dr. (82) 

If two fields are incident, e(r) and f(r), then we can define the maximum and 
minimum intensities as 



J \e(r) + f(r)\ 2 dr 

e(r)-f(r)\ 2 dr, (83) 
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and thus, 



Ima* = J (e(r) + f(r)) (e*(r) + /») dr 

= J (\e(r)\ 2 + \.f(r)\ 2 + .f(r)e*(r)+e(r)r(r)) dr 

= h + I f + J (f(r)e*(r) + e(r)/*(r)) dr. (84) 



Similarly, 



h + If-J (f(r)e*(r) + e(r)r 



Imin = I e + If- / (/(r)e*(r) + e(r)/* (r)) dr. (85) 



Therefore, the visibility of the interferometer is: 

v= J(f(r)e*(r) + e(r)r(r))dr _ 
h + If 

Of course, the relative phase between e(r) and /(r) was given by the definitions of 
the maximum and minimum intensities, and therefore /(r)e*(r) is real. Assuming 
that the total incident intensities of the two fields are equal, I e = If = I, the 
visibility of the interferometer is simplified to: 

v= ff(r)e*(r)dr _ (g?) 

This expression allows us to derive a number of relations, 

J J \e(x)f(y)\ 2 dxdy = J |e(x)| 2 dx J \f(y)\ 2 dy = I 2 (88) 

J J e*(x)f(x)e(y)f*(y)dxdy = J e*(x)f(x)dx J e(y) f (y) dy = I 2 V 2 (89) 

J j\e(x)f(y)-f(x)e(y)\ 2 dxdy = J J (e(x)f(y) - f(x)e(y)) 

x(e*(x)f*(y)-r(x)e*(y)) dxdy 
(\e(x)\ 2 \f(y)\ 2 + \f(x)\ 2 \e(y)\ 2 



-e*(x)f(x)e(y)f*(y) 
-e(x)f*(x)e*(y)f(y)) dxdy 
= 2I 2 - 2I 2 V 2 

= 2I 2 (1 - V 2 ) , (90) 

that will be useful below. 

The spatial modes £ appearing in the density matrix of ion b can now be ex- 
pressed in terms of the measureable parameters I and V. As above, we will assume 
the incident intensities are equal. 

P(\a\ 2 + \t3\ 2 (l-V 2 ) a P *V 2 \ 

Pb ~8\ /3a*V 2 |/3| 2 + M 2 (1 - V 2 ) J ■ iy > 
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Of course, the factor of I 2 /8 is part of the normalization factor neglected throughout 
this derivation; it is thereby discarded, and the condition tr(p b ) = 1 is used to 
determine the proper normalization. After doing so, we finally end up with the 
density matrix for ion b as: 



:: - ' ' |2 (1 - V 2 ) a(3*V 2 



Pb ~2-V 2 \ (3a*V 2 \(3\ 2 + \a\ 2 (l-V 2 )) ■ (92) 

The above enables calculation of the expected reduction in fidelity for any telcported 
state as a function of the visibility of the interferometer. For any given input state, 
the ideal density matrix is simply: 



/ \a\ 2 ad* \ 

Pideal=[y a ^^). (93) 



Therefore, the fidelity T = tr (pbpideai) is: 
T = tr 



1 f \a\ 2 + \(3\ 2 (1 - V 2 ) a/3*V 2 \ f \a\ 2 a(3* 

2-V 2 { (3a*V 2 \l3\ 2 +\a\ 2 (l-V 2 )) \/3a* \p\ 2 

= (M 4 + l/?l 4 + 2M 2 |/?| 2 (i - V 2 ) + 2\a\ 2 \P\ 2 V 2 ) 

= 2 ^7 2 (H 2 + I/3| 2 ) 
1 

~ 2-V 2 ' 



(94) 



We see that if the visibility of the interferometer is perfect, V = 1, then the fidelity 
is 1; whereas if the spatial mode-overlap is nonexistent, V — 0, then the fidelity 
drops to 1/2 (as expected for a totally mixed state). 

In the experiment, we measure the visibility of the interferometer by coupling 
laser light into the single-mode fibers used for transferring the spontaneously emit- 
ted photons from the atom to the beamsplitter. We find a visibility V > 0.98. By 
the above derivation, we can thereby estimate that the spatial mode-mismatch at 
the beamsplitter reduces the fidelity of the teleportation protocol by at most about 
4%. 



6.3.2. Imperfect state detection 

We calculate the expected degradation in fidelity due to imperfect state detection 
of the atomic qubit using the formalism of the measurement operators, as above. 
Let ej be the error in the measurement of the qubit state of ion j = a, b. The 
measurement operators may then be defined as: 

M 0j = yr^7|0)(0| + ^eJ|l)(l| 

^y = >/^|0><0| + v^^|l><l|- (95) 

Note that these satisfy the completeness relation Ij — ^2 m -M m j. 
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In this calculation we assume perfect interference of the photons at the beam- 
splitter, so that the state of the two ions after a coincident detection of the photons 



and rotation of ion a by R v {tt/2), yields the state (Eq. 64 1 



W)ions = ^|0)a (a|l> 6 + /3|0> 6 ) + ~|l) a (a\l) b - p\0) b ) 

= ~|0) (a\l) + (3\0)) + ^|1> («|1) " /?|0)) , (96) 

where in the second line we have again invoked the shorter notation of the prior 
section and suppressed the subscripts. 

The explicit derivation for measuring the state |0) a is presented below; the 
derivation for |1) Q is analogous, and after the conditional microwave rotation on 
ion 6, yields exactly the same result. As in the previous section, the measurement 
is completed by taking the partial trace over the measurement operator and the 
quantum state, so that following measurment of ion a the density matrix for b is 
given by: 



Pb = tr a [MoaMoal^ionaitp] 

= [Vn^Woi) (a\l) b + \0) b ) + V^UM) (o|l>6 - P\0) b )] 

® [Vl^a~\0)a (a* 6 (l| + (3* b (0\) + v^|l)a (<X* b {l\ ~ P* b( 

= (1 - e a ) (M 2 |l)(l| + /3a*|0)<l| + af3*\l}(0\ + |/3| 2 |0)(0|) 

+e a (|«| 2 |1)<1| - (3a*\0)(l\ - a/T |1)<0| + |/3| 2 |0)(0|) 
= |/5| 2 |0)<0| + (1 - 2e a )Pa*\0){l\ + (1 - 2e Q )a/3*|l)(0| + |a| 2 |l)(l|. (97) 

After measurement of ion a, a microwave pulse conditioned upon the measurement 
is applied to ion b. In the case of measuring |0) a , the rotation R x {tt) is applied to 
ion b. The density matrix is then: 



a\ 2 (l-2e a )ap* 
(1 - 2e a )l3a 



Notice that here the fidelity of the teleported state depends critically upon the state 
amplitudes. As can be seen in Eq. |98| if cither a or [3 is zero, then the imperfect 
detection of ion a plays no role in the final density matrix of ion b. Intuitively, this 
is correct because the influence of imperfect measurement of ion a is propagated 
by the conditional rotation on ion b. In the case that a or /3 is zero, then the two 
possible rotations R x {^) and R y (ir) perform the same action, and the measurement 
error of ion a is inconsequential. 

The fidelity of the operation with imperfect measurements can be evaluated 
by reconstructing the density matrix following measurement of ion b in the same 
fashion as is done with the actual experimental data. We define the probability of 
measuring a particular state \tp) as P\^)- Since all the measurements occur via state 
dependent fluorescence, measurement in bases other than the z-basis are completed 
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by performing a microwave rotation prior to detection. The six possible probabilities 
are then given by: 
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= tr 
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These probabilities allow us to calculate the Stokes parameters: 
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These parameters are the coefficients of the Pauli matrices in a simple analytical 
formula for the reconstruction of the density matrix^^ 

1 3 

Pb,recon = X ^ Sj&j. (101) 
j=0 

The reconstructed density matrix contains the effect of imperfect measurement 
on ion b. We can then calculate the fidelity of the teleported state, taking into 
account imperfect detection, by T — tr(pb,reconPideal)- Taking the state detection 
fidelities e a — 0.985 and e b = 0.975, we calculate the fidelity of the states { |0), |1), 
(|0) + \1})/V2, (|0) - |1))/ A (|0) + i|l))/A (|0) - i\l))/V2 } to be { 0.975, 0.975, 
0.961, 0.961, 0.961, 0.961 }, yielding an average reduction in fidelity of about 3.5%. 

Given that the fidelity of teleporting the two states |0) and |1) depends only 
on the measurement imperfection on ion 6, it would have been beneficial to have 
reversed the roles of the two ions in this experiment. Doing so could have improved 
the average fidelity by 0.003%. However, this change is within the uncertainty of 
the results presented here; the improvement would only be noticed in a much larger 
sample of events. 

6.4. Discussion 

As in the original teleportation proposal^ the successful implementation of our 
teleportation protocol requires the transmission of two classical bits of information: 
one to announce the success of the heralded quantum gate and another to determine 
the proper final rotation to recover the teleported state at ion b. While these classical 
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bits do not convey any information about the quantum states of either ion a or 6, 
in the absence of this classical information ion b is left in a mixed state (Eq. 65 1 , 
and the protocol fails. The required classical communication also ensures that no 
information is transferred faster than the speed of light. 

Despite the apparent differences, the heralded teleportation protocol reviewed 
her d 109 4s directly analogous to the original proposal from Ref.|7l In this comparison, 
consider ion b and photon b to be the entangled pair shared between Alice and Bob, 
as shown in Fig. [18] Thus, the three qubits used for the teleportation protocol are 
ion a, photon b, and ion b. Photon a is used in conjunction with the beamsplitter 
(BS) to perform an inefficient Bell-state measurement on ion a and photon b. 




Fig. 18. Comparison of the heralded teleportation protocoP^H and the original proposal from 
Ref.|3 The entangled qubits shared between Alice (system in the red, dotted box) and Bob (system 
in the blue, dashed box) are photon b and ion b. The required Bell state measurement between 
ion a (containing the information to teleport) and photon b is then performed (inefficiently) by 
generating photon a and interfering it with photon b at the beamsplitter (BS). The polarization 
filters (PBSs) allow only photons with the same polarization to be measured. The coincident 
detection of the two photons by the photomultiplier tubes (PMTs) heralds the success of the 
partial Bell state measurement. Measurement of ion a and a conditional single qubit rotation on 
ion b then completes the teleportation protocol. 



Although our implementation is intrinsically probabilistic because the four two- 
photon Bell states are not all deterministically distinguishablepIS the teleportation 
protocol succeeds without postselection due to the two-photon coincident detection 
that serves as a heralding eventP^The protocol presented here has the advantage 
of establishing the quantum channel between the (atomic) quantum memories using 
photons and entanglement swapping, allowing the atoms to be separated by a large 
distance from the outset. Ultimately, the teleportation scheme demonstrated here 
has the potential to form the elementary constituent of a quantum repeater capable 
of networking quantum memories over vast distances, and may be an essential 
protocol for the realization of scalable quantum computation. 
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7. Outlook 

The high fidelities obtained in the teleportation and quantum gate experiments is 
evidence of the excellent coherence properties of the photonic frequency qubit and 
the "clock" state atomic qubit. Together, these complimentary qubits may provide 
a robust system for practical quantum communication and quantum computation. 



7.1. Quantum gate with infrared photons 

In principle, any arbitrary distance can be bridged using a quantum repeater based 
on the entanglement protocols reviewed in prior sections. However, the number 
of nodes needed to efficiently implement the quantum repeater is approximately 
proportional to the inverse of the attenuation length of the photons! 3 ^ In order 
to establish quantum channels across long-distances, it may be more practical to 
use photons with wavelength in the infrared region of the spectrum (rather than 
the 369.5 nm photons used here), as these photons experience less attenuation in 
fiber. Moreover, access to additional optical frequencies may facilitate entanglement 
between disparate optically active systems, such as atoms and quantum dots. 

The rich atomic structure of the Yb + atom results in transitions across the op- 
tical spectrum. Two additional transitions that appear particularly amenable to 
the photon- mediated heralded gate described in Sec. [5] are the 935 nm 3 [3/2^/2 <-» 
2 Z) 3 /2 transition, and the 1.3 um. 2 P 3 / 2 2 Dz/2 transition. 

Spontaneously emitted photons at 935 nm whose frequency is entangled with the 



internal electronic states of the atom can be generated as illustrated in Fig. 19 'a). 
In this case, the atom is initialized in the 2 Si/ 2 ground state and excited using 
ultrafast pulses at 297.1 nm. Pulses at 297.1 nm could be generated by the third- 
harmonic generation of a mode-locked Ti:S laser operating at 891.4 nm. As long 
as the two atomic states are mapped to separate hyperfine manifolds in the 2 D 3 / 2 
level, then selection rules can be exploited to perform state dependent fluorescence 
detection of the atom as before; only population in the 2 D 3 / 2 \F — 1) manifold will 
be transferred to 2 Si/ 2 \F — 1) by the 935 nm light applied during detection. Given 
the 52 ms natural lifetime of the 2 D 3 / 2 level, it should be possible to obtain state 
detection fidelities > 98% P 



A heralded quantum gate can be performed using 935 nm photons (Fig. 19 a)) 
by first preparing the atom in a superposition of |0) and |1): 

\i/>)a = a\0)+p\l). (102) 

A 7r-polarized ultrafast pulse at 297.1 nm is used to coherently transfer the popu- 
lation from 2 5*1/2 to 3 [3/2] 1 / 2 - Due to the selection rules involved, this will transfer 
|0) to 3 [3/2] 1/2 |F=l,m F = 0) and |1) to 3 [3/2] 1/2 |F = 0, m F = 0). The 3 [3/2] 1/2 
level can then decay to 2 D 3 / 2 by spontaneously emitting a 935 nm photon. By col- 
lecting photons at 935 nm emitted perpendicular to the quantization axis, we can 
use polarization filters to distinguish w- and cr-polarized photons. A 7r-polarized 935 
nm photon is the result of a 3 [3/2] 1 / 2 |F = 1 , m p = 0) to 2 D 3 / 2 \F = 2,m F = 0) 
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Fig. 19. Quantum gate protocols utilizing infrared transitions in 171 Yb+. (a) Procedure for gener- 
ating photons at 935 nm with the frequency mode entangled with the 2 _D 3 /2 1^ 7 " = {1)2}, mp = 0} 
atomic state. As this protocol retains initial coherence in the atom, it is suitable for the im- 
plementation of a heralded quantum gate, (b) Method of generating 1.3 fim photons for the 
implementation of a heralded quantum gate. 

or 3 [3/2] 1 / 2 |-F 1 = 0,mp = 0) to 2 D 3 / 2 \F — l,mF = 0) transition, resulting in 
the frequency of the emitted photon being entangled with the internal state of the 
atom: 

\i>)a P = 0\F = 2, m F = 0)k) + a\F = l,m F = 0)\v b ). (103) 

Here, \F, mp) refer to states in the 2 D 3 / 2 level, and At/ = Vf, — v T = 3.07 GHz is the 
sum of the 2 D 3 / 2 and 3 [3/2] 1 / 2 hyperfine splittings. As the ion-photon entanglement 
process preserves the coherence initially present in the ion, the interference and 
detection of these photons can be used to implement a heralded quantum gatelMHH 
Detection of the atomic state can be accomplished by fluorescence detection at 369.5 
nm, as described above. In this case, detecting fluorescence at 369.5 nm indicates 
the atom was originally in the 2 D 3 / 2 \F — l,mp = 0) state, whereas the absence of 
fluorescence indicates the atom is in 2 D 3 / 2 \F — 2, mp = 0). 

The disadvantage of generating 935 nm photons is that the branching ratio 
of the 3 [3/2] 1 / 2 level between 2 Si/ 2 and 2 D 3 / 2 has been calculated to be about 
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55:1,^ decreasing the probability of generating a 935 nm photon and thereby re- 
ducing the overall success probability of the entanglment and gate protocols out- 
lined above. Due to this reduced success probability dark counts of the single- 
photon detectors could become a significant source of error. However, there is a 
way to "veto" the contribution of dark counts in this protocol. After a photon 
detection event, a microwave pulse at the 2 D 3 / 2 hyperfine splitting of 0.86 GHz 
can be used to transfer the 2 D 3 / 2 \F — l,mp = 0) atomic state to, for instance, 
the 2 D 3 / 2 \F = 2, trip = 1) state. If the fluorescence detection procedure outlined 
previously is now performed, no 369.5 nm photons should be detected]^] The popu- 
lation transfered from 2 D 3 / 2 \F — 1) to 2 D 3 / 2 \F — 2) can be returned by a second 
microwave pulse, and the remainder of the protocol completed. On the other hand, 
detection of 369.5 nm photons during this detection interval would indicate that 
the 935 nm "detection" was a false event, and should be discarded. 

The gate protocol outlined above could also be implemented at 1.3 fim, where 
attenuation in optical fiber is near a minimum (Fig. [l9^b)). In this case, the atom 
is initially prepared in 2 S\/ 2 , and then excited to 2 P 3 / 2 by an ultrafast pulse at 329 
nm. Decay from 2 P 3 / 2 to 2 D 3 / 2 results in the emission of a 1.3 /im photon. While 
the above protocol can be implemented in an analogous fashion at this wavelength, 
the branching ratio from 2 P 3 / 2 to 2 Si/ 2 versus 2 D 3 / 2 is about 475:lJ^Thus, while 
the wavelength is more amenable to long-distance transmission, the decrease in 
the protocol success probability is even more dramatic. In addition, 2 P 3 / 2 can also 
decay to 2 D 5 / 2 , which can subsequently decay to the long-lived 2 F 7 / 2 . Depopulating 
these additional metastable states would require additional optical frequencies and 
could limit the repetition rate of the experiment. 



7.2. Scalability 

The primary impediment to scaling the current setup (as well as the infrared proto- 
cols above) to more qubits is the success probability of the heralded quantum gate. 
As discussed in Sec. [5] there are several avenues that can be pursued to improve 
the success rate. Given that the effective solid angle of collection is about an order 



of magnitude smaller than the other contributing factors in Eq. 60 improvements 
here are bound to have the greatest impact. 

One proposal to increase the effective solid angle is to situate reflective or re- 
fractive optics near the trapped atomic ionl 115 l n £l A possible approach is to place 
the trapped ion at the focal point of a parabolic mirror, allowing the collection effi- 
ciency to approach unitylH- 7 l 118 llf properly segmented, a metallic mirror could also 
serve as the electrodes for the ion trap. Alternatively, the spontaneous emission into 
free space could be replaced by the induced emission into the small mode volume of 
a high finesse cavity, which can reach near unit efficiency) 119 ! 120 ! 121 ! 122 ! Even though 



°Here, the 14.7 GHz sideband usually used during Doppler cooling can be applied to the impinging 
369.5 nm light during this detection interval, to ensure population left in the |0) state also produces 
fluorescence. 
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the free spectral range of the cavity would have to be 14.7 GHz to simultaneously 
support both frequency modes in the gate protocol, choosing a near-concentric 
design could still result in a small mode volume and thus in a high emission prob- 
ability into a well-defined Gaussian mode. Of course, placement of optics near the 
ion will require careful assessment to ensure the added structures do not destabilize 
the quadrupole trap. As another option, the photonic qubit could be encoded in 
a temporal mode (time-bin qubit) , allowing more conventional cavity structures in 
the implementation of the heralded quantum gateP^ 

While improvements that increase the success probability of the gate operation 
can enhance scalability, even with a low success probability, this gate can still be 
efficiently scaled to more complex systems. The architecture for a quantum repeater 
could consist of photon-mediated operations between spatially separated nodes and 
local deterministic gates at each node. This implementation would also enable deter- 
ministic long-distance quantum teleportation and remote deterministic gates In 
this architecture, the photon-mediated operations between nodes can be attempted 
simultaneously The requirement for scaling to more nodes is therefore only that 
the coherence time of the atoms exceed the time needed to connect all the nodes 
of the quantum repeater. Given all connections can be attempted simultaneously, 
then the time needed to connect all nodes is approximately given by T success ln(iV), 
where T success is the average time needed to connect two nodes and N is the total 
number of nodes in the repeater (N large). With a measured coherence time of 2.5 
seconds for the 171 Yb + hyperfine qubitj^this means that to establish entanglement 
over a quantum repeater with 10 nodes, we require the average success rate to be 
greater than about 1 Hz(p] It should be possible to achieve this success rate with 
modest improvements in the photon collection efficiency and the repetition rate of 
the experiment. 

The stipulations on the success probability for generating large cluster states 
for scalable quantum computing arc more stringent. The time needed to construct 
an n qubit ID cluster state is approximately given by^ 

, , ^log 2 (n c + l) , v 

T(n)«t„ - + p-*— l°g 2 («-0, (104) 

\ * success J \ i success J 

where n c ss 4/P success is the critical number of qubits in a single ID cluster state 
chain that must be generated before multiple chains can be fused together. Clearly, 
even if we assume the repetition rate of the experiment can be improved to 100 ns, 
at our current success probability of 2.2 x 10~ 8 the time to generate a 100 qubit 
cluster state is prohibitively long (> 10 36 years). However, if the improvements 
suggested above are able to increase the gate success probability to 10%, then a 
100 qubit cluster state could be generated in less than 25 ms. 

p In this estimate, we have assumed only one photon-mediated connection is being attempted 
between each node. Multiplexing this operation further relaxes the required success rate. 
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7.3. Summary 

We have reviewed the recent experimental progress towards quantum networks 
with distant trapped atomic ions. We presented a detailed discussion of the recent 
implementation of a heralded quantum gate between remote ions based on the 
interference and detection of spontaneously emitted photons, and the employment 
of this gate to teleport a qubit between distant quantum memories. Ultimately, 
this quantum gate and teleportation protocol could be essential to realizing long- 
distance quantum communication and scalable quantum computation. 
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